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In this article we review the state of the art on the transport properties of quantum dot systems 
connected to superconducting and normal electrodes. The review is mainly focused on the theo- 
retical achievements although a summary of the most relevant experimental results is also given. 
A large part of the discussion is devoted to the single level Anderson type models generalized to 
include superconductivity in the leads, which already contains most of the interesting physical 
phenomena. Particular attention is paid to the competition between pairing and Kondo correla- 
tions, the emergence of 7r-junction behavior, the interplay of Andreev and resonant tunneling, and 
the important role of Andreev bound states which characterized the spectral properties of most of 
these systems. We give technical details on the several different analytical and numerical methods 
which have been developed for describing these properties. We further discuss the recent theo- 
retical efforts devoted to extend this analysis to more complex situations like multidot, multilevel 
or multiterminal configurations in which novel phenomena is expected to emerge. These include 
control of the localized spin states by a Josephson current and also the possibility of creating 
entangled electron pairs by means of non-local Andreev processes. 
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I. INTRODUCTION 

The field of electronic transport in nanoscale devices 
is experiencing a fast evolution driven both by advances 
in fabrication techniques and by the interest in potential 
applications like spintronics or quantum information pro- 
cessing. Within this context quantum dot (QD) systems 
are playing a central role. These devices have several dif- 
ferent physical realizations including semiconducting het- 
erostructures, small metallic particles, carbon nanotubes 
or other molecules connected to metallic electrodes. In 
spite of this variety a very attractive feature of these de- 
vices is that they can usually be described theoretically 
by simple "universal-like" models characterized by a few 
parameters. In addition to their potential applications, 
these systems provide a unique test-bed for analyzing the 
interplay of electronic correlations and transport proper- 
ties in nonequilibrium conditions. 

Electron transport in semiconducting QDs has been 
studied since the early 90's when pheno mena like 
Coulo mb blockade (CB) was first observed (Kastner[ 
1993). It soon became clear that QDs could allow to 



study the effect in transport properties of basic elec- 
tronic correlations phenomen a like the Kondo effect as 
suggested in early predictions (|Glazman and Raikh||T988l 
Lee an d Ng[ |l988| ). These predictions were first tested in 
metall ic nanoscale junctions cont aining magnetic impu- 
rities (Ralph and Burhman, 1994). However, a definitive 



breakthrough in the field came with the observation of 
this effect in semiconducting QDs by Goldhaber et al. 



Goldhaber-Gordon et al. 1998) and Cronenwett et al. 



Cronenwe tt et ar| "l998). A great advantage of these 



devices is to offer the possibility of controlling the rele- 
vant parameters, thus allowing a more direct comparison 
with the theoretical predictions. Since then the effect of 
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Kondo correlations in electronic transport has been ob- 
served in several physical realizations of QDs based on 
carbon nanotubes (CNTs) (Nygard et al. J2QQQ ) and big 



molecules like fullerenes (Roch et al. 2009 



In parallel to these advances the study of supercon- 
ducting (SC) transport in nanoscale devices has also ex- 
perienced a great development. From a theoretical point 
of view, with the advent of mesoscopic physics, a more 
detailed understanding of superconducting transport was 
developed around the central concep t of coherent An - 
dreev reflection (AR) (Andreev, 1964; Beenakker, 1997). 
This concept has allowed to unify the description of su- 
perconducting transport in different types of structures 
like normal metal-superconductor (N-S), S-N-S junctions 
and superconducting quantum point contacts (SQPC). 
Due to the multiple AR (MAR) mechanism the spectral 
density of systems like S-N-S or SQPCs is characterized 
by the presence of the so-called Andreev bound states 
(ABS) inside the superconducting gap. These states are 
sensitive to the superconducting phase difference and are 
thus current-carrying states which usually give the dom- 
inating contribution to the Josephson effect. 

In recent years it has become feasible to produce hy- 
brid systems combining different physical realizations of 
QDs well co ntacted to superconduct ing electrodes (for a 
review see (Franceschi et al., 2010)). Superconducting 



transport through QDs provides the interesting possibil- 
ity to explore the interplay of the AR mechanism and 
typical QD phenomena like CB and Kondo effect. The 
central aim of this review article is to discuss the main 
advances which have taken place on this issue during the 
last years. 

A usual assumption in these studies is that a basic 
description of the main properties of these hybrid sys- 
tems can be provided by the Anderson model and its 
generalizations to include SC leads, orbital degeneracy, 
etc. The single level model applies when the dot level 
spacing 5e is larger than all other relevant energy scales. 
In the normal state the model allows to describe in a 
unified way CB and the Kondo effect both in and out 
of equilibrium conditions fKouwenhoven and Glazman" 
2001| ). With two superconducting leads interesting new 
physics already appear in the equilibrium case. Due to 
the Josephson effect, electron transport is possible with- 
out an applied bias voltage and the model describes the 
competition between Kondo effect and induced pairing 
within the dot. Figure [l] illustrates this competition: de- 
pending on the ratio between the Kondo temperature, 
Tx, and the superconducting order parameter. A, there 
is a phase transition between a Kondo dominated spin- 
singlet ground state to a degenerate magnetic ground 
state. This transition is accompanied by a reversal of 
the sign of the Josephson current. Thus in the magnetic 
case the S-QD-S system constitutes a realization of the 



so-called 7r-junction (Glazman and Matveev 1989 Spi- 



vak and Kivelson, 1964). A more detailed understanding 






Source QD Drain 



Source QD Drain 



FIG. 1 Schematic representation of a single spin-degenerate 
level QD connected to normal (left panel) and superconduct- 
ing (right panel) leads with a large charging energy Uc- In 
the normal case the local density of states (LDOS) in the 
dot exhibits the typical form corresponding to the Kondo 
regime with a narrow resonance at the Fermi energy and a 
a broad resonance (of width F) below it. In the supercon- 
ducting case the Kondo resonance (assumed to be narrower 
than the superconducting gap) disappears due to the compe- 
tition with the pairing correlations in the leads. Courtesy of 
C. Schonenberger. 



cently (|Rozhkov and Arovas, 1999 Vecino et al., 2003). 
The realization of a 7r-j unction in QD systems should 
distinguished from the similar phenomena in SFS junc- 
tions, where F denotes a ferromagnetic material (Gol- 



ubov et al.\ 2004). 



of this transition describing the appearance of intermedi- 
ate phases with metastable states was achieved more re- 



Another basic situation which has been extensively ex- 
plored (both theoretically and experimentally) is the N- 
QD-S case. This situation has been mainly analyzed 
in the linear transport regime in which it exhibits and 
interesting interplay between Kondo behavior and reso- 
nant Andreev reflection. In contrast to the S-QD-S case 
this system does not exhibit a quantum phase transi- 
tion but there is instead a crossover from a Kondo domi- 
nated regime for large Tk / A to a singlet superconducting 
regime in the opposite limit. 

A third paradigmatic situation which has been stud- 
ied is the voltage biased S-QD-S system. This situation 
constitutes a much more demanding task for the the- 
ory due to the need of describing properly the strong 
out-of-equilibrium distribution which is generated by the 
infinite series of multiple Andreev reflection (MAR) pro- 
cesses together with the effects of Coulomb interactions. 
The problem becomes simpler in two limiting cases: 1) 
when Coulomb interactions are small and treated in a 
mean field approximation thus allowing to analyze the 
interplay of MAR and resonant tunneling and 2) when 
the Coulomb energy is the larger energy scale in the prob- 
lem and the contribution of MAR processes are largely 
suppressed. 

More recent developments include the study of several 
QDs (connected either in series or in parallel) coupled to 
one or more SC electrodes. In these situations there is 
a competition between not only the Kondo and the SC 
correlations but also the possible magnetic coupling of 
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the spins localized within the dots. This could open the 
possibility to control the spin state of the dots system by 
means of the Josephson current. 

In addition there is a growing interest in analyzing 
transport in these hybrid structures in a multiterminal 
configuration. One of the aims of these studies is the de- 
tection and control of non-local Andreev processes, which 
offer the possibility of producing entangled electron pairs 



strict the analysis to a single spin-degenerate level is 



(Recher et al 2001). 



This review article is organized as follows: in Section 
iniwe introduce the basic theoretical models used to de- 
scribe the hybrid QD systems, which are largely based 
on the Anderson model and its generalizations. In this 
section we also give a brief summary of the application 
of nonequilibrium Green functions techniques for the cal- 
culation of the electronic properties within these type of 
models. Section [1111 is devoted to review the main re- 
sults for the S-QD-S systems in equilibrium, i.e. in the 
dc Josephson regime. We first discuss the case of a non- 
interacting resonant level which is useful to illustrate the 
emergence of ABSs and its contribution to the Josephson 
current. In the subsequent subsections we give account 
of the different theoretical methods which have been used 
to analyze the effect of interactions in the dc Josephson 
regime. A main issue which is discussed in this section 
are the phase diagrams describing the transition to the 
TT-state as a function of the model parameters. We also 
give a brief account of the existing experimental results 
for S-QD-S devices in this regime. The case of a QD 
coupled to both a normal and a superconducting leads 
(N-QD-S) is addressed in Section IV Most of the re- 



sults obtained for this systems correspond to the linear 
regime with different levels of approximation to include 
the Coulomb interactions. We also briefly mention ex- 
isting results for the non-linear regime and the few ex- 
periments which have been reported of this case up to 
date. Section [V| is devoted to the voltage biased S-QD- 
S system. We first give some technical details on the 
calculations for the non-interacting case in order to il- 
lustrate how to deal with the out of equilibrium MAR 
mechanism. We also comment in this section the few 
existing results including interactions in this regime and 
give an account of the related experiments. Finally, in 
Section [Vl] we discuss several different situations which 
go beyond the single-level two-terminal case discussed in 
the previous sections. These include: multidot systems 
connected either in parallel or in series, the multilevel sit- 
uation and setups in a multiterminal configuration. We 
conclude this article with a brief discussion of related 
issues not included in the present review and of topics 
which, in our view, deserve to be further analyzed in the 
near future. 



II. BASIC MODELS AND FORMALISM 

The minimal model for a QD coupled to metallic elec- 
trodes in the regime where 5e is sufficiently large to re- 



provided by the single level Anderson model (Anderson 
1961), with the Hamiltonian H 



QD 



here Hqd corresponds to the uncoupled dot given by 



Hqd = ^ eo4a^Oa + Uno^nQ^, 



(1) 



where Cq^ creates and electron with spin a on the dot 
level located at eo and U is the local Coulomb interac- 
tion for two electrons with opposite spin within the dot 
(nocr = Cq^cqct)' On the other hand, Hl^r describe the 
uncoupled left and right leads which can be either nor- 
mal or superconducting. In this last more general case, 
they are usually represented by a BCS Hamiltonian of 
the type 

ka k 

(2) 

where ^ creates an electron with spin a at the single- 
particle energy level ^/c,zy of the lead v = L^R (usually re- 
ferred to the lead chemical potential, i.e. (^k,v = ^k,v — IJ^v) 
and = | Aj,| exp (i^^^) is the (complex) superconduct- 
ing order parameter on lead v. Finally, Ht describes the 
coupling between the QD level to the leads and has the 
form 



kcT^v 



(3) 



In order to reduce the number of parameters it is usu- 
ally assumed that the normal density of states of the leads 
pu{oj) is a constant in the range of energies around the 
Fermi level of the order of the superconducting gap and 
that the k dependence of the hopping elements Vku — 
can be neglected within this range. The coupling to 
the leads is then characterized by a single parameter 
Tjy = 7Tpu\Vy\^^ which can be interpreted as the normal 
tunneling rate from the dot to the leads. 

Within the above model Aj, = would correspond to 
the normal state. For vanishing Ty the model is in the so- 
called atomic limit which is characterized by sharp peaks 
in the spectral density at eo and eo + U . This limit cor- 
responds to the Coulomb blockade regime in an actual 
QD where the conductance is strongly suppressed except 
at the charge degeneracy points. When the couplings 
to the leads increase (F^^ become larger than tempera- 
ture) virtual processes allow the charge and spin in the 
dot to fluctuate and a resonance at around the Fermi 
energy appears close to half-filling due to the Kondo ef- 
fect. This simple model thus already captures the most 
relevant Physics of ultrasmall QDs with well separated 
energy levels, like the crossover from the Coulomb block- 
ade to the Kondo regime as the temperature is lowered. 

The simplicity of this model has allowed to obtain ex- 
act res ults in the equilibrium case by m eans of the Bethe 
ansatz (Wiegmann and Tsvelick, 1983). The most basic 
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of these results is the expression for the Kondo tempera- 



ture (Hewson, 1993) 



exp 



7r|eo(eo ■ 



2UT 



(4) 



where T = Tl^Tr. 

This temperature characterizes the crossover from the 
so-called local moment regime for T ^ Tk to the regime 
where Kondo correlations between the localized spin 
within the QD and the spin of the electrons in the leads 
sets in. Although this physics is basically well under- 
stood since the 70 's for the case of magnetic impuri- 
ties in metals, its consequences for transport in artificial 
nanostructures has started to be developed much more 
recently specially driven by the advances in fabrication 
techniques. In this respect while the linear transport 
properties are well understood still open questions re- 
main regarding the non-equilibrium regime. 

In the superconducting case another energy scale, as- 
sociated with the superconducting gap, appears bringing 
additional complexity to the problem, whose description 
is in fact the scope of this review. Even in the equilib- 
rium situation the Anderson model with superconducting 
leads contains the non-trivial Physics associated to the 
Josephson effect. A relevant parameter is then provided 
by the superconducting phase difference (j) = (1>l — (j^R- 

In order to analyze the electronic and transport prop- 
erties of a general superconducting system in the pres- 
ence of interactions and in a non-equilibrium situation 
it is convenient to use Green function techniques. The 
Keldysh formalism provides the basic tools for this pur- 
pose. 

Due to the presence of superconducting correlations 
it is convenient to introduce the Nambu spinor field op- 
erators ^j,^^, with = ^c^^,Cj|^ where j = kv^{) 

denotes the u = L^R electrodes and the dot level respec- 
tively. The different terms in the model Hamiltonian of 
Eqs. ([T]), ([2| and (|3| can then be written as 



Hqd = 

= 

Ht — 



k 

E(C^^-^o + h.c.) 



(5) 



where = cqTs, noa = [^o + sign(cr)r3] ^o, hku = 
^ki^Ts ^ReAj^Ti +ImAi^r2, r^^o, 1,2,3 being the Pauh ma- 
trices defined in Nambu space. 

Starting from these spinor field operators, general- 
ized single-particle propagators can be defined along the 
Keldysh closed time loop as 



straightforward way most of the relevant quantities like 
the mean charge and the superconducting order param- 
eter within the dot as well as the mean current through 
it, which are given by 



no{t) 
ReAo(t) 



zTr(r3GoV(^,^))-l (7) 
[/Tr(riGoV(^,^)) 

I ^ Tr (ra [T4.G+-o(t, t) - V^kG^oit^ t)] ) . 



These expressions are formally exact but of little use 
unless the single-particle propagators are known. Fully 
analytical and exact results can only be obtained in the 
non-interacting case. In the presence of interactions nu- 
merical methods allow to obtain exact results in the equi- 
librium case. In a more general case one is bound to find 
reasonable approximations for these propagators valid for 
a restricted range of parameters. It is usually convenient 
to express these approximations in terms of a self-energy 
S which is related to the propagators by the usual Dyson 
equation 



c(o) 



c(0)' 



Goo — Gqq + Gqq SooGoo, 



(8) 



:(o) 



where G denotes the unperturbed propagators cor- 
responding to an appropriately defined non-interacting 
Hamiltonian Hq (the "symbol indicates matrix structure 
in Keldysh space) and where integration over internal 
times is implicitly assumed. Different approximations for 
the self-energy associated with electron-electron interac- 
tions are discussed in the forthcoming sections. The anal- 
ysis of the problem is greatly simplified in the stationary 
case where Fourier methods can be applied both in the 
equilibrium and in the non-equilibrium situation. We 
shall start discussing in the next section the simplest pos- 
sible case of an equilibrium situation. In this case further 
simplification arises from the possibility of expressing all 
Keldysh propagators in terms of the retarded-advanced 
propagators and the Fermi equilibrium distribution func- 
tion, npioj) = 1/ [1 + exp P{oJ — /i)], where /i is the chem- 
ical potential and /3 = l/ksT as 



G+-(a;)=nFM[G^(a;)-G^(cj)] . (9) 
Thus, for instance, the mean current can be written as 



(6) 



where a, /3 = ± denote the two branches in the Keldysh 
contour. These propagators allow to calculate in a 



(10) 
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III. EQUILIBRIUM PROPERTIES OF QUANTUM DOTS 
WITH SUPERCONDUCTING LEADS 

For a nanoscale system coupled to superconducting 
electrodes a finite current can flow even in the absence 
of an applied bias voltage due to the Josephson effect. 



We shall illustrate this effect starting from the non- 
interacting situation {U = 0) within the single level 
Anderson model introduced above. In this case the 
advanced-retarded Green function of the coupled dot can 
be expressed as 



r 



ra. 
l9l 



(11) 



Leadl 



Nanostructure 
-M 



Lead 2 




FIG. 2 Schematic representation of the physical mechanism 
responsible for the formation of ABs in a generic nanostruc- 
ture coupled to superconducting leads. Reprinted by permis- 
sion from Macmillan Publishers Ltd: Nature Physics (Pillet 



et al 2010), copyright (2010). 



where /^'^ 



\^l,r\N\^l,r? - ('^ ± i0+)2 and g^' 



R 



1 L,R 

— {uj ± zO+)/"'^/|AL,i?| are the dimensionless BGS green 
functions of the uncoupled leads. For simplicity we fo- 
cus below on the case where both leads are of the same 
material for which \Al\ = \^r\ = A. 

The spectral density associated with this model ex- 
hibits bound states within the superconducting gap (i.e 
I Co' I < A). Physically, the ABSs arise from virtual multi- 
ple Andreev reflection processes at the interface between 
the dot and each of the leads. In such processes and for 
energies inside the gap, the electrons (holes) incident to- 
wards the leads are reflected back as holes (electrons) , as 
illustrated in Fig. |2j The condition for the appearance of 
ABSs is that the accumulated phase in the closed trajec- 
tory be a multiple of 27r, which is equivalent to satisfying 
the equation 



D{uj) = [uj-eo-Tg{uj)][uj^eo-Tg{uj)] 



where we have used that the dimensionless BCS Green 
functions f become real for energies inside the super- 



condu cting gap. It can be shown ( |Rozhkov and Arovas 



2000 ) that this equation has two real roots inside the gap 
icjs, i.e. symmetrically located with respect to the Fermi 
energy. 

An essential property of the ABSs is that they cor- 
respond to current-carrying states. In fact, due to 
its dependence on the superconducting phase differ- 
ence they have associated a Josephson current is(^) = 
2e/fi{duOs/d(j)). The total Josephson current is obtained 
by adding the contribution of all states with finite occu- 
pation. Thus, at zero temperature only the lower ABS 
contributes and there is an additional contribution from 
the continuous spectrum < — A which we discuss be- 
low. 

The ABS equation ([l2| becomes particularly simple 
for an electron-hole and left-right symmetric case (i.e. 
eo = and = F^^) when it can be reduced to 



UJ ± A cos (j)/2 ■ 



0, 



(13) 



(12) 



which for F A has the simple solutions ujs — 
±A cos 0/2, where A is a reduced gap parameter which 
for |(^| < 1 is given by A = A [l - 2(A/F)2] . The ABSs 
for this case tend to the ones of a perfectly transmit- 
ting one channel superconducting contact ±A cos 0/2, 
the main qualitative difference being the reduced am- 
plitude of their dispersion detaching them from the gap 
edges at = 2n7r. This is illustrated in Fig. |3|a). 

In the absence of electron-hole symmetry (lc. eo ^ 0) 
a finite internal gap between the upper and lower ABSs 
appears as in the case of a non-perfect transmitting one 
channel contact. When A/F ^ the ABSs for this 

case are given by uJs = iA^^ 1 — r sin^ (0/2), where 

r = 1/(1 + (eo/F)^) is the normal transmission at the 
Fermi energy. Outside this limiting case the ABSs ex- 
hibit both the internal gap and the detachment from the 
continuum states at = 2n7r, as it is illustrated in Fig. 
ib). 

An interesting issue to comment is that the contribu- 
tion from the states in the continuous spectrum becomes 
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FIG. 3 Andreev bound states for a non-interacting S-QD- 
S with different F/A values: 1.0, 2.0 and 4.0. Left panel 
corresponds to the e-h symmetric case (eo = 0) and the right 
panel to a case with eo = 0.5 A. 



continuum states. The expression for the total Josephson 
current in the non-interacting case can be derived from 
Eq. (10), which yields the following compact form 



h 



Fi,Fi^sin(/) / dcjni? (ct;)Im 



£r -Cr 

JlJ R 

]jr 



(15) 



The current density in Eq. ( 15 ) contains both the con- 
tribution from the ABSs (region Ico'l < A) and from the 
continuous spectrum |cc;| > A. The behavior of the cur- 
rent density as a function of A/F is depicted in Fig. [4] 
As can be observed the contribution from the continuous 
spectrum has the opposite sign compared to the one aris- 
ing from the ABS. The inset shows that this contribution 
becomes negligible in the limit A/F ^ and reaches a 
maximum for A/F ~ 2. 

In the rest of this section we shall discuss the different 
theoretical approaches to include the effect of interac- 
tions in the dc Josephson effect through single level QD 
models. We also include a subsection on experimental 
results. 



A. Cotunneling approach 




FIG. 4 Current spectral density for the non-interacting S- 
QD-S system with eo and (j) — 1.5 for increasing values of 
F/A = 0.5 (red), 4 (green) and 16 (blue). A small broadening 
has been introduced to help visualizing the ABs contribution, 
which has been truncated for the sake of clarity. The inset 
shows the relative contribution of the continuous spectrum Ic 
compared to the total value Ij. 



negligible when A/F ^ 0. In this limit the zero temper- 
ature current-phase relation (CPR) is simply given by 



^s(0) 



eA 



r sm ( 



2% 



1-rsin^ (0/2) 



(14) 



which is the CPR of a one channel contact with transmis- 
sion r. For finite A/F there is a contribution from the 



From a theoretical point of view the simplest approach 
to account for the effect of interactions in the Josephson 
current is to perform a perturbative expansion to the 
lowest non-zero order in the tunnel Hamiltonian. This 
so-called cotunnelin g approach was first used by Glaz- 
man and Matveev (Glazman and Matveev, 1989[ ), who 
predicted the onset of the 7r-junction behavior by this 
method. More precisely they obtained for the U ^ oo 
limit 



h A 



sm ( 



(16) 



where A changes its value from 2 (eo > 0) to -1 (eo < 0), 
thus describing the transition to the 7r-phase, the func- 
tion F{x) having the form 



F{x) 



dtidt2 



(coshti + cosht2)(^ + coshti)(x + cosht2) 

(17) 

This approximation is clearly not valid for describ- 
ing the Kondo regime {Tk ^ A) which requires non- 
perturbative approaches like the ones discussed in fol- 
lowing subsections. 



B. Mean field and variational methods 

Another simple approximate methods to deal with in- 
teractions are those of a mean field type like the Hartree- 
Fock approximation (HFA) or the slave-boson mean field 
(SBMF). In spite of their simplicity these approximations 
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are able to capture important qualitative features due to 
interactions in certain limits. 

We start by analyzing the HFA. In the context of mag- 
netic impurities in superconductors this method was first 
applied by Shiba ( ,Shiba^ ^1973j), while for the analysis 
of the Josephson effect it was first considered in Ref. 
( Rozhkov and Arovas , 1999) and further analyzed in (Ve-| 
cino et al. , 2003 ) . Within this approximation electrons 



with a given spin "feel" a static potential due to the aver- 
age occupation of the opposite spin, which corresponds to 
a simple constant self-energy (Hqo)!! = U < noi > and 
(^00)22 = —U < no^ >. In principle within the same 
level of approximation there appears a non-diagonal self- 
energy taking into account the induced pairing within 
the dot due to proximity effect, which can be written as 
(^00)12 = U < cj^cj^ >. The effect of this non-diagonal 



contribution, whic h was not included in Ref. ( Rozhkov 
and Arovas |1999 ), was analyzed in Ref. (Vecino et al. 



2003). The determination of the self-energy in the HFA 
requires a self- consistent calculation by using Eqs. (|8| 
which cannot be performed analytically in general. 

The most significant result within the HFA is the ap- 
pearance of a broken symmetry state in which the dot ac- 
quires a finite magnetic moment (i.e. < uq^ >7^< ^oi >) 
for certain ranges of parameters. In this respect one 
should be cautious in principle as the HFA is known to 
predict also broken symmetry states for the same model 
with normal leads (Anderson] [1961 ), which are known 
to be spurious. However, for the superconducting case 
ground states with a finite magnetization do exist for cer- 
tain parameter range as commented in the introduction. 
As it is shown below the HFA gives a rather good esti- 
mate of the magnetic ground state energy in the regions 
where it is the most stable phase. 

The general properties of the ground state within the 
HFA are most conveniently displayed by a phase-diagram 
like the one in Fig. [5] In this diagram the notation "0", 
"0^", "tt^" and "tt" corresponds to the different ground 
state symmetries. Thus, "0" corresponds to the non- 
magnetic case for all values of (j) (the absolute energy 
minimum being located at = 0), while the "tt" denotes 
that the magnetic solution is the most stable for all (j) 
values (with the absolute minimum at ^ = tt). On the 
other hand, "0^" and "tt^" refer to intermediate situa- 
tions with mixed magnetic and non-magnetic solutions 
as a function of (/), the name indicating whether the ab- 
solute energy minimum corresponds to a non-magnetic 
or a magnetic solution. From Fig. |5] the broken sym- 
metry ground states are predicted to appear around the 
eo = —U/2 line, which corresponds to the half- filled case 
for sufficiently large [/, i.e. U > F, A. It is worth notic- 
ing that for normal leads this region corresponds to the 
deep Kondo regime, which anticipates an interesting in- 
terplay between both effects in the superconducting case 
beyond the HFA. 

Further insight into the HFA solution can be provided 



by a "toy" model introduced in Ref. (Vecino et al. , 2003 



(A similar model was analyzed in Ref. (Benjamin et al. 




FIG. 5 Phase diagram in the U, — ep plane for F/A 



tained using the HFA (Rozhkov and Arovas 1999)) 



1 ob- 

^ The 

phases are classified into 0, , tt and tt as explained in the 
text. Reprinted figure with permission from A.V. Rozhkov 
and D. Arovas, Physical Review Letters 82, 2788, 1999 



(Rozhkov and Arovas 1999). Copyright (1999) by the Amer- 
ican Physical Society. 



2007)). In this simplified model the finite magnetiza- 
tion which appears in the HFA is simulated by means 
of an exchange field parameter, E^x^ corresponding to 
the splitting of the diagonal dot levels for each spin, i.e. 
^Ocr = Co + crEex- The analysis of the Andreev states 
within this toy model is similar to the one given at the 
beginning of this section for the non-interacting case, and 
becomes particularly simple in the limit A <C F in which 
they adopt the analytical expression 

f^±Y - ^^^^ + + ^^(^^ + sin^ ± 2XS{(I)) 



\ AJ 



Z4 + 2(X2+^2)^1 



(18) 



where E = eo/2r, X = Eex/2T and 
S{(j)) is given by 



S{(l)) = \JZ^ cos2 (t)/2 + + sin^ ^/4. 

This expression clearly shows that the effect of the ex- 
change field is to break the spin degeneracy producing 
an splitting of the ABSs. Consequently for Eex 7^ one 
could in principle observe up to four ABSs in the spec- 
tral density. The evolution of these states with increasing 
Eex is shown in Fig. |6] together with the corresponding 
Josephson current. While for E^x < T the splitting is 
small and states corresponding to different spin orienta- 
tion do not cross, for increasing E^x the upper and lower 



8 



<1 



< 



0.0 

-IjO 
1.0 

0^ 

0.0 

-0.5 

■1.0 
1.0 

0.5 

0.0 

-0.5 





/ \ 




/I ® 




f 



become 



O.S0 




0.25 








0.00 


1 










-0,25 








'0.50 




0.50 








o.as 




0.00 


1 






-025 




-0-50 




0.50 








0.25 




0.00 


1 










-0,25 




-0.50 





IjO 

<t)/7C 



0,5 1,0 1,5 



FIG. 6 Evolution of the ABSs and Jose phson current in the 
"toy" model of Ref. ( jVecino et al]^003l for eo/F = -0.5 and 
increasing Eex/^ parameter: 0.25 (upper panels), 0.75 (mid- 
dle panels) and 1.5 (lower panels). Reprinted figure with per- 
mission from E. Vecino et aL, Physical Review B 68, 035105, 
2003 dVecino et aH|2003t . Copyright (2003) by the American 
Physical Society. 



states closer to the Fermi energy begin to cross yielding 
a current-phase relation of 0' or tt' character. Eventually 
for sufficiently large E^x these two states completely in- 
terchange position with a complete reversal of the sign 
of the current (7r-phase). It should be mentioned that 
although in this toy model the spin degeneracy is artifi- 
cially broken, it nevertheless qualitatively simulates the 
behavior of the exact solution in the magnetic phase. 

Another simple approach of a mean field type is 
provided by the slave boson mean field approximation 
(SBMFA). This method was introduced by Coleman for 



the normal Anderson model (Coleman 1984 Hewson 



1993| ). It is based in the introduction of auxiliary bo- 



son fields ^0 5^0 which act as projectors onto the empty 
impurity state. At the same time fermion creation and 
annihilation operators /q^, /ocr are introduced for describ- 
ing the singly occupied states. In order to get rid of the 
doubly occupied states in the U ^ oo these operators 
should satisfy the completeness relation 



(19) 



Hqd 



y^go/ocr/Ocr 
a 

E(^^4.,.^S/o<. + h-c.) 



(20) 



So far this transformation is exact in the U ^ oo 
limit. Specific diagrammatic methods to obtain the im- 
purity self-energy in this slave boson formulation have 
been developed (Bickers, 1987). Within this formulation 



the simplest solution is provided by the mean field ap- 
proximation in which the boson operator is treated as a 
c-number. In fact the dot Hamiltonian reduces in this 
case to 

H^D = E 'ofLfo. + A (^bo\' + E fLfo. - , 

(21) 

where A is a Lagrange multiplier associated to the con- 
straint (19). The problem becomes equivalent to a non- 



interacting impurity model with renormalized parame- 
ters eo = eo + A and = l^opr^^. Self-consistency is 
achieved by minimizing the system free energy. 

Strictly speaking the mean field approximation is only 
valid in the N ^ oo limit where N is the level degener- 
acy of the Anderson model (for the single level case N=2). 
However, the mean field approximation yields a reason- 
ably good description of qua ntities like the Kondo tem- 
perature in the normal case (Hewson, 1993). When ap- 



plied to the Anderson model with superconducting elec- 
trodes the SBMFA is only valid in the regime Tk ^ A 
as it is not able to describe the transition to the 7r-phase 
when Tk ^ A. In the regime Tk A the ABSs 
as described by the SBMFA corresponds to the non- 
interacting case with renormalized parameters eo and 
f. In this way the ABSs within the SBMFA in this 

regime would be given by ujs{^) = ±^1 — f sin^ (/)/2 with 

f = 4T|^/(eQ + 4T|^). In principle, the self-consistent ef- 
fective parameters in the superconducting state can differ 
from those in the normal state. However, in the limit 
Tk a in which the approximation is supposed to 
work this difference can be neglected. The SBMFA in 
the U ^ oo limit has only been applied for the case of 
superconducting lea ds in a few references: Avishai et al. 
( Avishai et a/. J2003) for analyzing the dc current with an 



applied bias voltage, and in Ref. (Yeyati et al 2003) for 
studying the dynamics of Andreev states in the Kondo 
regime. Both works correspond to the non-equilibrium 
situation which will be discussed in Section [Vl 

For a proper description of the phase-diagram within a 
mean- field slave boson approach a finite- version of the 
method, like the one introduced by Kot liar and Ruck- 
enstein dKotliar and Ruckensteinl 1986), is necessary. 



In terms of these operators the terms Hqjj and Ht 



Within this method the number of auxiliary boson fields 
is extended up to four, denoted by e^Pa- and (i, which 
project into the empty, singly occupied (with either spin 
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orientation) and doubly occupied dot states respectively. 
These operators must satisfy the constraints 



E 



ptPa 



- e^e ■ 



pIpo 



(22) 



For recovering the non-interacting limit it is necessary 
to introduce also an auxiliary operator Za- = {1 — d'^ — 
^2^-i/2^g^^ +p^(i)(l — — p|)~"^/^, in terms of which 
the model Hamiltonian becomes 



<J kv',a 



(23) 



where A and Act are the Lagrange multipliers associated 
with the constraints (22). Again, in a mean field approx- 



imation, the auxiliary fields are treated as c-numbers to 
be determined self-consistently. 

The type of phase-diagram that is obtained within the 



finite-U SBMFA will be analyzed in subsection III.D[ for 
the zero band- width limit which allows a comparison with 
exact diagonalizations. As it is shown in that subsection 
the finite-U SBMFA tends to underestimate the stability 
of the TT-phase in contrast with the HFA, which typically 
overestimates it. 

Another relatively simple approach is provided by the 
use a variational wave-functio n. This a pproach was used 
by Rozhkov and Arovas (Rozhkov and Arovas^ 200 ) ex- 
tending previous works ( Varma and Yafet , 1976 ) in which 



variational wave-functions were proposed for analyzing 
the normal Kondo problem. In their work Rozhkov and 
Arovas propose different many-body variational states in 
the U ^ oo limit for the singlet and the doublet states, 
looking for the their relative stability. They find a transi- 
tion between both ground states for A/Tk ^ 2 and also 
predict the appearance of the intermediate phases 0' and 
7t' in addition to the pure and tt ones. 
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FIG. 7 Feynmann diagrams for the second-order self-energy 
in the Anderson model with superconducting electrodes. 
Reprinted figure with permission from E. Vecino et a l, Phys- 
ical Review B 68, 035105, 2003 f Vecino et a/.||2003|. Copy- 



right (2003) by the American Physical Society. 



C. Diagrammatic approaches 

1. Perturbation theory in the Coulomb interaction 

A natural extension over the HFA is provided by ap- 
plying diagrammatic perturbation theory in the Coulomb 
parameter U. Already at the level of second order one 
can obtain an approximation for the self-energy which 
is able to capture part of the interplay between Kondo 
effect and pairing interactions. This approximation has 
been applied both for the S-QD-S case in equilibrium 



( Matsumoto , 20011 Vecino et al 2003 ) , as well as for the 



N-QD-S case (Cuevas et a/., SOOTfYamada et al 2010). 



The diagrams contributing to the second-order self- 



energy in the superconducting Anderson model are de- 
picted in Fig. [t] The first diagrams (denoted as 11(a)) 
is equivalent to the one appearing in the normal case, 
describing interaction of an electron with an electron- 
hole pair with opposite spin. The other diagrams include 
anomalous superconducting propagators and are there- 
fore characteristic of the superconducting state. The 
presence of these propagators gives several effects: the 
appearance of non-diagonal elements of the self-energy 
in Nambu space, and the presence of diagrams like 11(b) 
and 21(b) in Fig. [t] which corresponds to a double- 
exchange process. Finally, diagram 21(a) describes the 
interaction of a Cooper pair with fiuctuations in the pair- 
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ing amplitude within the dot. 

Formahy, these diagrams can be computed from the 



fuh Green functions of the non-interacting case by means 
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The evaluation of these expressions for a general range 
of parameters requires a significant numerical effort. An 
efficient algorithm can be implemented to evaluate these 
expressions based on Fast Fourier transformations, as dis- 



cussed in (Schweitzer and Czycholl, 1990). 



In the limit A/F <C 1 Kondo correlations dominate 
over pairing ones. The results of the second-order self- 
energy approach for the half-filled case capture the main 
features of the onset of Kondo correlations in the spectral 
density whe n U > T. Th i s is il lustrated in Fig. [S] taken 
from Ref. ( jVecino et al\ 2003), which shows its evolu- 
tion for increasing U /F. The spectral density is similar to 
the one found in the normal state except for the super- 
imposed features inside the superconducting gap. The 
overall shape evolves from the single Lorentzian broad 
resonance for < F to the three peaked structure char- 
acteristic of the Kondo regime when U > V. In this 
regime the width of the central Kondo peak is set by the 
scale Tx, which in the present approximation is given by 
Tk F/(1 — Qfo), where 



behavior of the non-interacting case but with a nar- 
rower dispersion set by a;s(0) A [l — ([//[/q)^], where 
(f/o/r)2 = (r/A)7rV(27r + 2). 

On the other hand, when Tk ~ A a transition to the 
TT-phase is expected. Within the second-order self-energy 
approach the transition can be identified by allowing 
for a breaking of the spin-symmetry in the initial non- 
interacting problem and searching for self-consistency. 
Rather than imposing the consistency co ndition of the 
HFA, i.e. eocr = eo + ^ < ^Oa > in Ref. (IVecino et al 



2003 ) it was imposed that the effective dot level for each 



spin-orientation be determined by the charge-consistency 
condition, i.e. < noa >=< Uq^ >, where < Uq^ > is the 
dot charge corresponding to the broken-symmetry non- 
interacting Hamiltonian. Such a procedure was shown to 
eliminate the unstable behavior of perturbation theory 



when developed from the HFA (Yeyati et al 1993) 



2. NCA approximation 



(0) 



TT 



(24) 



thus coinciding with the perturbative result in the normal 
state (see Ref. (Yamada and Yoshida , 1975| )). Although 
this perturbative approach fails to yield the exponential 
behavior of Tk for large f7/F, it provides a reliable de- 
scription of the spectral density fo r moderate values of 
this parameter (Ferrer et a/.[ |1986 ). 



The second-order self-energy allows also to analyze 
the renormalization of the ABSs due to the presence of 
Coulomb interactions. For values of U/T < 10 the renor- 
malized ABSs maintain approximately the ~ cos<p/2 



Within the diagrammatic approximations one can in- 
clude the so-called non-crossing approximation (NCA). 
In this case an infinite order resumation of the perturba- 
tion theory is performed starting from the U ^ oo slave 
boson rep resentation of the Anderson Hamiltonian (B ick-] 



ers 



1987| ). To the lowest order in 1/N (where N = 2 is 
the spin degeneracy) the family of diagrams in this re- 
sumation is represented in Fig. |9] The dashed lines cor- 
respond to the fermion propagators and the wavy lines 
to the slave bosons. In the normal case the NCA in- 
clude only the first two diagrams in the Dyson equation 
for the fermion and boson propagators. The extension 



to the superconducting case was proposed in Ref. ( Clerk 
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FIG. 10 Josephson current and subgap LDOS for the equilib- 
rium S-QD-S model within the NCA for three values of Tk / A 
and three different temperatures. Reprinted figure with per- 
mission from G. Sellier et aL, Physical Review B 72, 174502, 
2005 (Se llier et a/.[|2005t . Copyright (2005) by the American 
Physical Society. 



FIG. 8 Evolution of the DOS in the S-QD-S system in equi- 
librium within the second-order self-energy approximation for 
and electron- hole symmetric case with A/F = 0.1. The U/F 
parameter takes the values 2.5 (a), 5 (b) and 10 (c). Reprinted 
figure with permission from E. Vecino et ai, Physical Review 
B 68, 035105, 2003 ( | Vecino et am2003J . Copyright (2003) 
by the American Physical Society. 



n(iv) = 




FIG. 9 Fermion (top) and boson (bottom) self-energy dia- 
grams in the NCA approximation extended to the supercon- 
ducting case. Reprinted figure with permission from G. Sell- 
ier et al, Physical Review B 72, 174502, 2005 ( [Sellier et aJ] 
[2Q05p . Copyright (2005) by the American Physical Society. 



Clerk et. al. (iClerk and Ambegaokar, 2000) is also for- 
mally exact to order 1 /N and it is thus expected to yield 
reasonable results even in the presence of MAR processes. 

Further analysis of the N CA applied to the S-QD-S 
system was provided in Ref. (Sellier et al. 20051). Their 



results for the Josephson current and LDOS in the super- 
conducting gap region are summarized in Fig. 10 The 



fact that the calculations are performed for temperatures 
which are a quite large fraction of A yields very broad 
resonances for the subgap states. As can be observed 
in Fig. [Toj only one broad resonance can be clearly re- 
solved within the gap. These results are in contrast to 
what is obtained using exact numerical methods as will 



be discussed in Sect. III.D.2 In this approximation the 
transition to the 7r-phase appears as a smooth crossover 
which can be associated to the crossing of this resonance 
through the Fermi energy. 



3. Real time diagrammatic approach 



et al 2000 ) and corresponds to including the anomalous 



propagators for describing multiple Andreev reflection 
processes (last two diagrams in the fermion and boson 
self-energies represented in Fig. [9|. 

In the normal case the NCA theory has been shown 
to yield reliable results for tem peratures down below Tk 
(Bickers 1987 Hewson, 1993|) in spite of certain patholo- 
gies like its failure to fulfill the Friedel sum-rule. The 
self-consistent extension for the superconducting case by 



Another technique which has been applied to the study 
of quantum dots coupled to superconducting leads is 
the real time diagrammatic appr oach first introduced by 
Konig et al. (Konig et al, 1996| ) for the normal Ander- 
son model. The main idea of this technique is to inte- 
grate out the fermionic degrees of freedom of the elec- 
trodes leading to a reduced description of the density 
matrix projected on the Hilbert space of the isolated dot 
states. In the superconducting case this reduced density 
matrix also depends on the number of Cooper pairs in 
the leads relative to some chosen reference. The aim of 
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the technique is to determine the time evolution of this 
reduced density matrix in the Keldysh contour thus al- 
lowing to consider both equilibrium and non-equilibrium 
situations. In Ref. (Governale et al. , '2008) the method 



A further step in the idea truncating the Hilbert space 
is performed in the so called zero band- width limit (Af- 



has been applied to the equilibrium S-QD- S case obtain- 
ing results in agreement with those of Ref. ( Glazman and 
Matveev, 1989") in the cotunneling limit. The method has 



been mainly applied to analyze the properties of quan- 
tum dots connected to both normal and superconducting 
leads in multiterminal configurations out of equilibrium, 
an issue which will be commented in Sect. I VII 



D. Diagonalization by numerical methods 

Within this subsection we will review methods which 
attempt a direct diagonalization of the superconducting 
Anderson model, either by truncating the initial Hilbert 
space using physical arguments valid for certain param- 
eter region or by using the Numerical Renormalization 
group (NRG) method. 



1. Exact diagonalization for the large A limit 

An exact diagonalization of the model is possible in the 
limit A ^ oo. In this limiting case the Hilbert space of 
the problem is automatically reduced to states spanned 
by the different electronic configuration of the dot lev- 
els. The effect of the superconducting leads appears as a 
pairing term between the electrons within the dot. The 
effective Hamiltionian for the truncated Hilbert space be- 



comes (T. Meng and Simon, 2009 Tanaka et al] 2007b 



Vecino et al] 2003) 



Jjeff ^ 2r COS ^/2 (co^coi + 4i4t) + ^0 

cr 

^Uno^no^. (25) 

The eigenvalues of this reduced Hamiltonian can be 
determined straightforwardly by noting the decoupling of 
subspaces with even and odd number of electrons. The 
ground state for the even case (corresponding to total 
spin /S = 0) is a linear combination of the empty and 
doubly occupied dot state with an energy 



Eois=o){<l>) = eo + U/2 - ^{eo + U/2y + AT^ cos^ <l>/2. 

(26) 

On the other hand, the odd number subspace sim- 
ply corresponds to a single uncoupled spin with energy 
^0(5=1/2) = ^0- The transition to the magnetic ground 
state thus occurs for £^o(5=i/2) = ^o(5=o)(</^)- the 
simpler electron-hole symmetric case (eo = —U/2) this 
condition reduces to 2rcos^/2 = U/2 and thus the full 
TT state appears for F < [//4. This simple model already 
gives a rough qualitative account of the — tt quantum 
phase transition. 



fleck et aL\ [2QQQ1 [Bergeret et~al] |2QQ7| [Vecino et al 



2003). In this approximation the superconducting leads 



are represented by a single localized level (which formally 
corresponds to the limit of vanishing width of the leads 
spectral density). This approximation is justified when 
the superconducting gap is large compared to the other 
energy scales in the problem, and thus can be consid- 
ered as a refinement with respect to the previous ap- 
proach. The corresponding Hamiltonian can be written 
8iS H = Hd^ Ht ^ Hl ^ Hr, with 

= ^ ^vc\^Cycj + (A^c^^c^,; + h.C.) 

cr 

Ht = ^ (KcLco. + h.C.) , (27) 

where u = L^R denotes the left-right sites describing the 
leads in this approximation. 

Although the total number of particles is not a good 
quantum number, their parity is conserved as in the pre- 
vious case. This allows to reduce the initial 64 states in 
the Hilbert space to a subspace of 20 states for even par- 
ity with a total spin z-component 5'^ = and 15 states 
for odd parity with Sz = ±1/2. These values of Sz are 
the ones corresponding to the ground state in each sub- 
space with total spin 6* = and S = 1/2. In addition 
to providing a qualitative description of the phase dia- 
gram of the full model, this simplified calculation can 
furthermore be useful as a test for comparing different 
approximation methods. 

The phase diagram obtained within this approxima- 



tion was dis cussed in Refs. ( Bergeret et al. 2007 Vecino 
et al] |2003D and is shown in Fig. |ll|for F = F = A, 
with Vl = Vr = V. As can be observed, the overall dia- 
gram is very similar to the one shown before for the HFA 
(Fig. [5| exhibiting the four phases 0,0^7r' and tt in the 
same sequence. For a more direct comparison it is nec- 
essary to perform the HFA of the ZBW model, which as 
an exactly solvable model also provides a stringent test 
of the approximation. The lower broken line in Fig. 11 
indicates the boundary of the 7r-phase within the HFA. 
It can be noticed that the HFA overestimates the sta- 
bility of this magnetic phase. On the other hand, it is 
also possible to test the finite-U SBMF approximation in 
this ZBW model. The corresponding boundary for the 
TT-phase is indicated by the upper broken line in Fig. [TT 
In opposition to the HFA this approximation overesti- 
mates the stability of the phase, the exact boundary 
therefore lying in between the two different mean-field 
approximations. It is interesting to point out that a sim- 
ilar difference between both approximations is also found 
for the full model (see inset in Fig. 11). One would then 
expect that the exact boundary for the full model lays in 
between these two. 

The ZBW model can also be used to test appr oximate 
methods beyond the mean field ones. In Ref. (Vecino 
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FIG. 11 Phase diagram of the S-QD-S system in the ZBW 
model for the leads obtained by exact diagonalization for 
A = Vl = Vr, taken as unit of energy. The dashed lines cor- 
respond to the boundary between the and tt regions within 
the HFA (lower hne) and the finite-U SBMFA (upper line). 
The inset show the corresponding results for these bound- 
aries within the full model. Reprinted figure with permission 
fromF.S. Bergeret et a /., Physical Review B 76, 174510, 2007 
(Bergeret et al.\ |2007|. Copyright (2007) by the American 



Physical Society. 



et ( 



2003 ) this was done for the second-order self-energy 



approximation. As it is illustrated in Fig. 12 for the sym- 
metric case (eo = —U/2) the second-order self-energy ap- 
proach matches quite well the exact ground state energy 
for U/ A values up to 2.5. It is interesting to note that 
for larger U values even when the HFA already yields 
a full TT-state, the second-order approximation predicts a 
mixed ground state in agreement with the exact solution. 
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FIG. 12 Ground state energy for the 5' = and 5' = 1/2 
states of the S-QD-S system in the ZBW model for the leads. 
Full lines correspond to the exact results, dotted lines to the 
HFA and dashed lines to the second-order self-energy approx- 
imation. Reprinted figure with permission from E. Vecino et 
al, Physical Review B 68, 035105, 2003 (I Vecino et a/. 112003). 



Copyright (2003) by the American Physical Society. 



Yoshioka and Ohashi (Yoshioka and Ohashi, 2000) and 



implemented by several authors to analyze the S-QD-S 



2004 



model with a finite phase difference ( Choi et al. 
Lim and Choil [20081 |Qguri et al \ |2QQ4| |Tanaka" et al 
2007b). For the left-right symmetric case (i.e. Al = 
Aji and F l = ^r) the sequence of Hamiltonians can be 



written as ( Choi et al 2004 ) 



/i,cr 

h. c. (^8) 



2. Numerical Renormalization Group (NRG) 

The NRG method is based on the ideas of Wilson on 
logarithmic discretization for magnetic impurity prob- 



lems (Wilson 1975 ) and was first applied to the Anderson 
model by Krishna- murthy et al. ( Krishna- murthy et al. 



[1980). The idea behind the method is to discretize the 
energy levels in the leads on a logarithmic grid of energies 
A~^ (with A > 1 and 1 < n < N ^ oo) with exponen- 
tially high resolution on the low-energy excitations. This 
discretization allows then to map the impurity model into 
a linear "tight-binding" chain with hopping matrix ele- 
ments decaying as A~^/^ with increasing site index n. 
The sequence of Hamiltonians which is constructed by 
adding a new site in the chain is then diagonalized iter- 
atively. As the number of states grows exponentially an 
adequate truncation scheme is required. 

The NRG scheme has been first generalized to the case 
of an Anderson impurity in a superconducting host by 



where the initial Hamiltonian is given by 



Ho 



/A 



^Ocr/M,0,cr + h. C. 



fi—e,o a 



(29) 



The fermion operators f^^N.a correspond to an effec- 
tive tight-binding chain resulting from the logarithmic 
discretization and the canonical transformation into the 
even-odd linear combination of original left-right states 
in the leads and 




A =A 



cos^/4 , Vo 



(30) 



2r 

-XA/ r7^sin^/4, 
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FIG. 13 Phase diagram of the S-QD-S system in the Tr^cq 
plane for fixed U and different values of Fl/Fh obtained using 
the NRG method. Reprinted figure with permission from Y 
Tanaka et a/., New Journal of Physics 9, 115, 2007 (Tanaka 
[7t a/.||2007b| ). Copyright (2007) by lOP Publishing Ltd. 



with X = 2/(1 + 1/A) and D being an energy cut-off in 
the leads spectral density. The original Hamiltonian is 
recovered in the limit H/D = lim^oo Hjsf / {x^^^~'^^^'^)• 
Th.e NRG method was applied to analyzin g the Joseph- 
son c urrent in a S-QD-S system in Ref. ( |Choi 



et al 



2004), this work confirming the predicted — tt quan- 
tum phase transition at A ~ Tk for the electron-hole 
symmetric case. It should be mentioned that more re- 



cent calculations ( [Karrasch et aL\ |2QQ8, ) using NRG ob- 
tain Josephson currents whic h are approximate ly a fac- 
tor 2 larger than the ones of (Choi et ar||2QQ4l). On the 



other hand, Oguri et al. ( |Oguri et al] |2QQ4| ) used NRG 
to analyze this model in the case of \Al\ In 
this case the model can be exactly mapped into a sin- 
gle channel model consisting on the right lead coupled to 
the Anderson impurity with a local pairing A^^ = Tlc'^^^ , 
thus allowing a simpler implementation of the NRG algo- 
rithm. Further work for the Tl ^ Tr case a lthough with 
Al = Aji by Tanaka et al. (Tanaka et al , |2QQ7b[ ) con- 
firmed the presence of intermediate 0' — tt phases even 
in the left-right asymmetric case. A characteristic phase 



diagram obtained for this case is shown in Fig. 13 



In addition to the ground state properties, NRG meth- 
ods have been applied in an attempt to clarify the struc- 



ture of the subgap ABSs. In Ref. (Lim and Choi 2008) 



the spectral density inside the gap obtained from the 
NRG algorithm was analyzed, showing that a pair of 
ABSs located symmetrically respect to the Fermi energy 
is present in the U oo case. This is in contrast to 



the NC A results discussed previously (shown in Fig. 10 ) 
where a single broad reson ance appears. Similar conclu- 
sions are obtained in Ref. (|Bauer et a/., 2Q07[ ) although 
for the single lead case and for finite U. A word of caution 
should be said regarding the analysis of the ABSs in this 



o 
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FIG. 14 Current phase-relations for the S-QD-S system ob- 
tained using the fRG approach truncated at the HF level for 
different values of U/T and A/F. For comparison the results 
obtained using the NRG method are also plotted (indicated 
by the filled dots). Reprinted figure with permission from C. 
Karrash et al, Physical Review B 77, 024517, 2008 (Karrasch 



et al] 2008). Copyright (2008) by the American Physical So- 



ciety. 



last work in which the relation Ti22{^) = is as- 

sumed in their Eq. (8) for the states inside the gap. This 
relation would not be strictly valid for the doublet ground 
state when choosing a given spin orientation. In this case 
the quasi-particle excitation energies would become spin 
dependent and the electron-hole symmetry would be bro- 
ken. This would allow in principle to have up to 4 ABSs 
inside the gap as predicted both by the Hartree-Fock ap- 
proximation and in the exact A ^ oo limit. Of course, 
in the tt phase the spin is not frozen but is fluctuating. 
In this sense the above relation between the self-energy 
components would be recovered when averaging over the 
Sz = 1/2 and Sz = —1/2 states. We believe in any case 
that a more detailed analysis of the ABSs using the NRG 
method is still lacking. 



E. Functional renormalization group 

The functional renormalization group (£RG) method is 
based on the application of an RG analysis to the dia- 
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grammatic expansions in terms of electronic Green func- 
tions. This is an approximate method whose accuracy 
depends on the initial diagrams used in the evaluation 
of the electron self-energies. The starting point is the 
introduction of an energy cut-off A into the Matsubara 
non-interacting Green-functions 

G^^^{iuj) = G{\uj\-A)G\iuj) 

Using these propagators the n— particle vertex functions 
acquires a A dependence. The flow equations are deter- 
mined differentiating these vertex functions with respect 
to A which are then solved iteratively for increasing A. In 
Ref. ( Karrasch et al. 2008 ) the method has been applied 



to the S-QD-S system employing a truncation scheme 
which keeps only diagrams corresponding the the static 
Hartree-Fock approximation. The A-dependent Green 
function used in Ref. ( Karrasch et al. , 2008 ) was of the 
form, 



luj — eo 



A{iuj) - icD + eo + 



(31) 



where iuj = iuj{l + X]/^ ^^^(ia;)) and A{iuj) = 
r^/(zcj)e*'^^, g and / being the dimensionless BCS 
Matsubara Green functions of the uncoupled leads. 
Within this approximation the flow equations lead to 
energy- independent self-energies, corresponding to an ef- 
fective non-interacting model with renormalized param- 
eters. It is important to notice that this approximation 
exhibits also the limitation already pointed out in the 
previous section as it imposes electron-hole symmetry 
which is not satisfied in the magnetic phase. Neverthe- 
less the approximation allows to identify a transition to 
a phase with inversion of the Josephson current which is 
driven by an " overscreening" of the induced pairing de- 



termined by Sa- Fig. 14 shows the comparison of fRG 
results with those obtained with the NRG method. The 
agreement is rather good for large A/F but it becomes 
poorer in the 7r-phase. We believe that the agreement 
could be improved allowing for a broken symmetry state 
within the same fRG approach. 



F. Quantum Monte-Carlo 

The Josephson current in the S-QD-S system has also 
been analyzed using Quantum Monte Carlo (QMC) simu- 
lations by Siano and Egger (Siano and Egger, 2004). The 
method used was the Hirsh-Fisher algorithm adapted 
to this particular problem. They consider the deep 
Kondo regime U/T 1 and eo/F <C —1 and show that 
the results for the Josephson current exhibit a univer- 
sal dependence with Tk provided that U/T > 5. They 
identify the transition between the different phases at 
A/Tk 0.51,0.875 and 1.105 for - 0^ 0^ - tt' and 
7t' — 7T respectively (Siano and Egger |2005a). Being a fi- 



relations do not exhibit sharp discontinuities in the in- 
termediate phases. This smooth behavior was criticized 



in Ref. (Choi et al 2005a) pointing out that the QMC 



results did not match the NRG ones of Ref. ( Choi et al. 



2004) at finite temperatures, which was attributed by 



Siano and Egger in their reply ( Siano and Egger 2005b ) 



to a possible limited accuracy of the NRG calculation of 
Ref. (Choi et a/. , 2004). More recent NRG calculations of 



Ref. ( Karrasch et al. , 2008| ) give a good agreement with 
QMC results at finite temperature for (j) values between 
7r/2 and tt, whereas QMC underestimates the current in 



the ra nge — 7r/2. It is claimed in Ref. ( [Karrasch et al.^ 
[2008 ) that the origin of the discrepancy lies in the fact 



that the first excited state in this phase range is smaller 
or of the order of the temperat ure va lues used in the 
calculations of (Siano and Egger, 2004). 



The QMC method has more recently been applied to 
analyze the spectra l prop erties of this model in Ref. 
(Luitz and Assaad 2010). The authors employ the 



so-called weak-coupling continuous-time version of the 
method which is based on a perturbative expansion 
around the U = limit. They show that the results for 
the spectral densities are in qualitative good agreement 
with the ones obtained in the zero band-width approx- 
imation introduced in Ref. (Vecino et a/., 2003), which 



was discussed before in this section. 



G. Experimental results 

Several physical realizations of the S-QD-S system have 
been obtained in the last few years by means of contact- 
ing carbon nanotubes (CNT), C60 molecules or semicon- 
ducting nanowires with su perconducting electrode s (for 
a recent review see Ref. ( [Franceschi et al. 2010[ )). In 
view of the existence of this review on the experiments 
in this section we give only a brief summary of the main 
findings and its relation to the theoretical work. 

CNTs have provided so far the most promising setups 
for a direct test of the theoretical predictions concerning 
the Josephson effect through a QD. The first experiments 
detecting a supercurrent through a CNT-QD strongly 
coupled to the leads (i.e. F ^ A, C/) were performed by 
Jarillo-Herrero et al ( Jarillo-Herrero et al. |2006 ). These 
experiments were basically performed in the resonant- 
tunneling regime with a single- level spacing ^ F. The 
results indicated a strong correlation between the criti- 
cal current Ic and the normal conductance Gn- However, 
the product IcRn deviated from a constant value due to 
the effect of the electromagnetic environment suppressing 
the critical current more strongly in off-resonance condi- 
tions due to phase fluctuations. 

In a subsequent work by the same group the first exper- 
imental evidence of 7r-j unction behavior in a S-QD-S sys- 
tem was obtained using semiconducting InAs nanowires 



with Al leads in a SQUID configuration ^van Dam et al. 
2006). However, the observed features in this experi- 



nite temperature calculation the resulting current-phase ment could not be explained completely on the basis of 
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FIG. 16 Color scale plots of the local density of states in a 



CNT QD coupled to SC leads as reported in Ref. (Fillet et al 



2010 ) as a function of gate voltage and phase difference. The 
figure also shows the results obtained from a phenomenologi- 
cal model, as discussed in the text. Reprinted by permission 



from Macmillan Fublishers Ltd: Nature Fhysics ( [Fillet et al. 
20101, copyright (2010). 



FIG. 15 Experimental setup used in Ref. ( [Cleuziou et al.\ 
2006 ) for analyzing the — tt transition in CNT QDs (upper 
panel). The middle panel shows the reversal of the oscil- 
latory pattern of the current as a function of the magnetic 
flux through the SQUID across the transition. The lower 
panel shows the correlation of the vr-phase region with the 
Kondo ridges in the normal state (indicated by the dashed 
line). Reprinted by permission from Ma cmillan Fublishers 
Ltd: Nature Fhysics ( [Cleuziou erarj|2006| ), copyright (2006). 



a single-level model but rather a multilevel description 
was necessary. In particular the authors showed that 
in this case the 7r-j unction behavior is not necessarily 
linked to the parity of the number of the electrons in the 
dot (this will be further discussed in Sect. VI). On the 
other hand, 7r-j unction behavior have also been demon- 



strated in Ref. ([ Cleuziou efal] |2006[ ) using CNT-QDs 
in a SQUID geometry. The corresponding experimental 
setup is depicted in Fig. 15 In this SQUID configuration 



the transition to the 7r-phase is directly demonstrated in 
the measured current-phase relation as a function of one 
of the applied gate voltages (see Fig. 15). Remarkably, 
this experiment allows to correlate the appearance of the 
TT-j unction behavior with the presence of Kondo correla- 
tions in the normal state. The dotted lines in the lower 
panel of Fig. [15] indicate the Kondo ridge appearing in 
the normal state (which is recovered by applying a mag- 
netic field). 

Evidence of a — tt transition has also been found in 
CNT-QD systems by analyzing the cur rent- volt age char- 
acteristic (Jorgensen et al, 2007). In this work it was 
shown that the evolution of the zero-bias conductance 
can be correlated to the behavior of the critical current 
when transversing the — tt boundary. As it corresponds 
to a non-equilibrium situation this analysis will be dis- 
cussed later in Sect. [Vl A similar observation holds for 



other experimental works ( Eichler et al] 2QQ9| Grove- 
[Rasmu ssen et a/.[ |2QQ7 ; Vecino et aL\ 20041 which will be 
discussed in that section. 
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Finally, it is worth pointing out recent experimental 
developments which have allowed to directly measure the 
Andreev bound states spectrum of a CNT-QD coupled 
to superconduct ing leads in a SQUID configuration ( |Pil-| 



let et al. 2010). In this experiment a weakly coupled 



lead was deposited at the center of the CNT to allow 
for tunneling spectroscopy measurements. In this way 
both the phase and gate voltage variation of the An- 
dreev bound states were determined. The results could 
be fitted satisfactorily using a simplified phenomenolog- 
ical model corresponding to the superconducting Ander- 
son model within a mean field approximation, in which 
an exchange field is included to represent the magnetic 
ground state. An example of the comparison between 
theory and experiment is given in Fig. [T6j This analysis 
showed that a double dot model was in general neces- 
sary to fit the experimental results. Recent experimental 



work on graphene QDs coupled to SC leads (Dirks et al 



[2011 ) provided also evidence on the crossing of ABs as a 
function of the gate potential which is consistent with a 
magnetic ground state. 



IV. QUANTUM DOTS WITH NORMAL AND 
SUPERCONDUCTING LEADS 

A single-level quantum dot coupled to a normal and a 
superconducting lead provides a basic model system to 
study electron transport in the presence of Coulomb and 
pairing interactions. Compared to the S-QD-S situation, 
this case has a simpler response in non-equilibrium con- 
ditions due to the absence of the ac Josephson effect. For 
this reason this system has been widely analyzed theo- 
retically. 

As in any N-S junction the low bias transport proper- 
ties are dominated by Andreev processes. This mecha- 
nism is in general highly modified by resonant tunneling 
through the localized levels in the dot. In addition, charg- 
ing effects can strongly suppress the Andreev refiection 
in certain parameter ranges. Furthermore there is also 
an interesting interplay between Kondo and pairing cor- 
relations as in the case of the S-QD-S system discussed 
above. 

As an illustration of the general formalism we derive 
here the linear transport properties of the non-interacting 
model. One can straightforwardly write the dot retarded 
Green function for this case from expression (11) by set- 
ting Al = 0. Then, from the expression of the current 
in terms of Keldysh Green functions and using the cor- 
responding Dyson equation one can write 



II 



E 



dcuTi 



(32) 



where g'^j^ ' ^ are the Keldysh Green functions of the 
uncoupled normal lead. By further using GJq"' ^ = 



+- 



and taking the wide band 
approximation for the uncoupled leads one can obtain 
the following expression for the temperature dependence 



linear c onductance ( Cuevas et al. 
mondi[|1999^ 



2001 Schwab and Rai- 



G 



16e^ 
h 



TnTs / dujlm^Gl^Gl^) 



driF 
duo 



(33) 



which, at zero temperature reduces to the simple expres- 
sion arising from the contribution of pure Andreev refiec- 
tion processes 



G = 



4e2 



S 



(34) 



As shown in Ref. ( [Beenakker 1992 ) this expression for 
the non-interacting case is equivalent to the formula G = 
(4eV/i)TV(2 - t)\ with T = ATnTs/{4 + (F^ + Tg?) 
being the normal transmission through the dot at the 
Fermi energy. In contrast to the case of a NS quan- 
tum point contact with essentially energy independent 
transmission, in the dot case the Andreev processes be- 
come resonant at eo = reaching the maximum value 
G = Ae^/h. 



A. Effect of interactions (linear regime) 

One of the first attempts to describe the effect of 
Coulomb interactions in the NDQS system was presented 
by Fazio and Raimondi (Fazio and Raimondi [1998 ) us- 
ing the equation of motion technique truncated to the 
second order in the tunneling to the leads. They derived 
expressions for the mean current using the K eldysh for- 
malism and extending the so-called Ng ansatz (Ng 1993) 



to the superconducting case. The claim of an extended 
temperature range for the zero bias anomaly due to the 
Kondo resonance was subsequently corrected in (Fazio 



and Raimondi 



1999j). 



The problem was addressed in Ref. ( |Kang , 1998) by 



assuming that the relation G = (4e^//i)r^/(2 — r)^ of 
the non-interacting case still holds by substituting r by 
the normal transmission of the interacting case. Within 
this assumption that work suggested an increase of the 
conductance in the Kondo regime by a factor of two with 
respect to the normal case. As shown by subsequent 
works which we discuss below, this enhancement is not 
always possible, the general case is rather the opposite. 

The conductance of th e interacting N-QD-S system 
was also analyzed in Ref. (Schwab and Raimondi, 1999) 
within the infinite-U slave-boson mean field approxima- 
tion. This approximation reduces the problem into an 
effective Fermi liquid description with renormalized pa- 
rameters Tn,s and eo. Within this approximation both 
Tn^s are renormalized equally and therefore the condi- 
tion for the maximum conductance is reached for the 
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symmetric case as in the normal state. As the authors 
acknowledge this result is valid only in the deep Kondo 
regime A <C T^, otherwise residual interactions would 
renormalize the left and right tunneling rates differently, 
in particular Ts coupling the dot with the superconduc- 
tor would be significantly suppressed by interactions. 
The problem was subsequently addressed by Clerk et 
^000 ) using the extension of the NCA 



al. (Clerk et 



to the superconducting case already mentioned in Sect. 
|III[ They analyze three different models: the N-QD-S 
Anderson model and a single channel magnetic and a 
two-channel non-magnetic contact between normal and 
superconducting electrodes. We comment here only the 
results for the first model. They find an overall decrease 
of the quasiparticle spectral density at the Fermi energy 
together with the appearance of additional Kondo peaks 
at ±A. As a consequence, their conclusion was that there 
is no enhancement of the linear conductance due to An- 
dreev processes in contrast to the claim of previous works. 
Diagrammatic techniques for the finite-U N-QD-S An- 



derson model were applied in Ref. (Cuevas et a/., 2001). 
Within this approach the linear conductance can be ex- 
pressed as the one corresponding to the non-interacting 
model with asymmetrically renormalized parameters 



G 



4e2 



4p2 f 2 



(^1 + r^ + f|)^ 



(35) 



where = - 1112 (0) and eo = eo - Sn (0) , 11^,^, being 
the dot self-energy elements in Nambu space. Although 
this result is valid in general within a diagrammatic anal- 
ysis, concrete results were obtained in this work by means 
of an interpolated second-order approach. From Eq. (35) 



it was found that even when the starting bare parameters 
correspond to the symmetric case Fat = F^, interactions 
would tend to reduce the conductance by inducing an 
asymmet ry, i.e. leading to Ts ^ Ijy, as was suggested 
in Ref. (Schwab and Raimondil 19991). However, this 



equation also predicts the possibility that an adequate 
tuning of the bare coupling parameters could yield an 
enhancement of the conductance up to the unitary limit 
{Ae^ /h in the NS case) which would not correspond to 
the maximum conductance in the normal case. 

The approximation of Ref. (Cuevas et al. 2001) is 



based on the evaluation of the second order diagrams, 
which due to the proximity induced pairing in the dot 
are formally the same as those of Fig. [7| which were dis- 
cussed for the S-QD-S case. For the non-symmetric case 
the authors used an interpolative ansatz which recovers 
the correct behavior in the F/[/ ^ (atomic) limit. The 
obtained behavior of the conductance as a function of 
U /F for the symmetric and non-symmetric cases is illus- 



trated Fig. 17 



As can be observed, in the symmetric case the con- 
ductance drops steadily from the unitary limit as U /V 
increases, the scale of this decay being set by the param- 
eter F/A. On the other hand the right panel of Fig. [17^ 
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FIG. 17 Linear conductance for the N-QD-S system in the 
secon d-order self-energy approximation of Ref. ( Cuevas et al. 
2001). (a) symmetric case as a function of U /V for different 



values of F/A = 0.125, 0.25, 0.5, 1.0, 2.0, 4.0 and 8.0, from top 
to bottom, (b) Same as (a) for asymmetric coupling to the 
leads Fl / Fh and different values of U/Vl (within the text 
we have set Vl =Tn and F^ = F5). Reprinted figure with 
permission from J.C. Cuevas et a l, Physical Review B 63, 
094515, 2001 ( Cuevas et a/.[[200T] ). Copyright (2001) by the 
American Physical Society. 



illustrates that the unitary limit can be restored by an 
adequate tuning of the ratio Fat/F^. 

The behavior of the linear conductance in the N-QD- 
S was also analyzed in (Domanski and Donabidowicz| 



2008 Krawiec and Wysokinski| |2004p using the EOM 

While in 



techni que with different decoupl i ng sc hemes . 
Ref. ( Krawiec and Wysokinski[ 2004) it was obtained 
that the conductance due to Andreev processes was com- 
pletely suppressed in the Kondo regime, an improved ap- 
proximation for the EOM decouplin g procedure in Ref. 
(Domanski and Donabidowicz , 2008) showed that there 



is in fact a finite zero bias anomaly in the Andreev con- 
ductance although it is in general much smaller than the 
one in the normal case, its value depending on the ratio 
Tat/F^. This behavior is in qualitative agreement with 
the results of the diagrammatic approach discussed be- 
fore. However, the unitary limit is not reached within 
this approach. 

More recently the linear conductance of the N-QD-S 



model has be en studied using the NRG method ( Tanaka 
|2QQ7a). In this work the limit A ^ 00 was taken 



et 



from the start, which allows to map the problem into the 
case of a QD with a local pairing amplitude A^^ = F^ 
coupled to a single normal electrode (the A 00 limit 
for the S-QD-S case was discussed in Sect. III). The NRG 



algorithm can be considerably simplified by this assump- 
tion because a simple Bogoliubov transformation allows 
to get rid of the local pairing term leading to a prob- 
lem which is formally equivalent to a normal Anderson 
model, which implies Fermi liquid behavior. As in the di- 
agrammatic approach discussed before, the main effect of 
the interactions is to renormalize the couplings Tn^s and 



the dot level position. Figure 18 illustrates the evolution 
of the renormalized parameters Fat, 5 as a function of U 
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FIG. 18 Renormalized couplings to the leads Fat, 5 (upper 
panel); Induced pairing amplitude (middle panel) and linear 
conductance (lower panel) as a function of U/Tn for the N- 
QD-S system obtained using NRG method in the A ^ 00 
limit. Reprinted figure with permission from Y. Tanaka et 
al, Journal of the Physical Society of Japan 76, 074701, 2007 
( Tanaka et arj[2007a| ). Copyright (2007) by the Physical So- 



ciety of Japan. 



for the half- filled case cq = —U/2 with initial parameters 
F^ = SFat. As can be observed in the upper panel of 
Fig. ~ 




FIG. 19 Color-scale map of the linear conductance in the 
N-QD-S system in th e Ts /U-cq/U plane obtained in Ref. 
( [Tanaka et al] |2007a| using the NRG method as in Fig. [I8 
Reprinted figure with permission from Y. Tanaka et a /.. Jour- 
nal of the Phy sical Society of Japan 76, 074701, 2007 ( [Tanaka 
et al] |2007a|. Copyright (2007) by the Physical Society of 



Japan. 



half-filled case obtained from these NRG calculations is 
illustrated in Fig. 
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In this color-scale map it can be 
clearly observed that the unitary limit is reached mainly 
along the curve {{eo/U + 1/2)^ + {Ts/U)Y^^ = 1/2, 
which corresponds to the single-doublet transition in the 
Fat ^ limit. When Ts/U < 0.5 (corresponding to 
the doublet state in the F^v ^ limit within the dashed 
curve in Fig. [l9|, the conductance as a function of eo ex- 
hibits a double peaked structure, whereas for Ts/U > 0.5 
only a single peak is found which can be correlated to the 
superconducting singlet ground state of the system in the 
Fat ^ limit. 



B. Non-linear regime 

The non-linear regime in the N-QD-S system has re- 
ceived so far much less attention and there are still as- 
pects, specially those related to the Kondo effect which 
are not sufficiently understood. This regime has been 
analyzed using the EOM technique employing different 



decoupling schemes in Refs. (Fazio and Raimondi 



is to reduce the effective coupling to the superconduc- 
tor, F^, while the effective coupling to the normal lead 
Ftv remains almost constant up to the region U ^ lOT^ 
when the Kondo effect is significant. Eventually the sys- 
tem reaches the condition F^r — and the conductance 
increases up to the unitary limit, as shown in the lower 
inset of Fig. 



Krawiec and Wysokinski[ |2 004^ 'S un et al] |200l| and 



1998 



TS) the main effect of increasing the interaction U ([Domanski and Donabidowicz, 2008D. The approxima- 



tion used in Ref. (Fazio and Raimondi 1998) has been 



already described in the context of the linear regime. On 
the other h and, the approximation used in Ref. (Sun 



et al, 2001) consisted in truncating the EOM equations 
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This behavior is in good agreement 
with the prediction of the diagrammatic theory of Ref. 
(Cuevas et al. 2001). On the other hand, the induced 



pairing amplitude (middle panel in Fig. 18), exhibits the 



expected monotonous decrease for increasing intradot re- 
pulsion. 

The behavior of the linear conductance outside the 



at the level of the two particle Green functions by sub- 
stituting the leads fermionic operators by their average 
values. Within this decoupling the authors find the ap- 
pearance of Kondo features in the dot spectral density. 
In addition to the usual features at = ieF for the 
regime A > F^ > F^v, they also find excess Kondo like 
features at = ±(2eo -\- U — eV). These features have 
been explained as arising from co-tunneling processes in- 
volving Andreev tunneling from the QD-S interface and 
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FIG. 20 Non- linear conductance for the N-QD-S system ob- 
taine d using the EOM method in the d ecoupling scheme of 
Ref . ( [Domanski and Donabidowicz 2008 ) in the A ^ oo limit 
as a function of cV/Vn and Fs/FAr. Reprinted figure with 
permission from T. Domanski and A. Donabidowicz, Physi- 
cal Review B 78, 073105, 2008 (iDomanski and Donabidowiczl 



2008) . Copyright (2008) by the American Physical Society. 



normal tunneling from the N-QD interface. 



In Refs. (Domanski and Donabidowicz 2008 Krawiec 



and Wysokinski , 2004 ) already commented for the linear 
regime, the finite voltage case was also analyzed. The 
non-linear conductance obtained for the A ^ oo limit 



in Ref. (Domanski and Donabidowicz 2008) both as a 



function of the bias voltage and the asymmetry in the 
coupling parameters F^/Fat is shown in Fig. 20 The 
parameters of this case correspond to the Kondo regime 
of the normal state eo = — I.SFat and U = lOF^r. This 
figure exhibits the expected features like the qualitative 
evolution of the zero bias anomaly already discussed in 
the previous section and the splitting of the dot reso- 
nances due to the proximity effect. Additional peaks at 
eV ^ ±U can be observed corresponding to the popula- 
tion of the higher charge states. 

The non-linear case has been analyzed more recently 



in Ref. (Yamada et al. 2010) by extending the interpola- 



tive self-energy approach of Ref. (|Cuevas et al\ |2001|) to 
the non-equilibrium case. The authors consider the case 
A = Fat and eo = —U/2 for different values of F^ and 
U . Their main findings are illustrated in Fig. 21 corre- 
sponding to the non-linear conductance for U = 2OF7V. 
They observe a Kondo peak which is displaced from zero 
bias, whose height increases with increasing F^ while its 
position is only weakly modified. When reducing F^ a 
second peak develops which shift progressively towards 
the gap edge cV/Tn = 1- 

We sho uld a lso mention the work of Ref. (Koert 
2010| ) in which the non-equilibrium transport 



ing et 



through a N-QD-S system is studied within an effective 
cotunneling model. Within this approach the self-energy 
is calculated to leading order in the cotunneling ampli- 
tude from which the nonlinear cotunneling conductance 
can be obtained. By neglecting charge fluctuations in 
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FIG. 21 Non- linear conductance for the N-QD-S system us- 
ing the inte rpolative self-energy approach of Ref. { Yamada] 
et al.\ 2010). The results correspond to the symmetric case 
with Fat = A and different values of Ts and U. Reprinted 
figure with permission from Y. Yamada et aL, Journal of the 



Physical Society of Japan 79, 043705, 2010 (Yamada et al 



2010). Copyright (2010) by the Physical Society of Japan. 



the dot two different regimes are found corresponding to 
the case of even and odd number of electrons. For the 
even case the system becomes equivalent to an effective 
S/N junction with the subgap transport due to Andreev 
reflection processes. On the other hand, for the odd case 
they find that the net spin within the dot leads to the 
appearance of subgap resonances giving rise to a peak- 
dip structure in the differential conductance. The typi- 
cal conductance curves that are found for both cases are 
shown in Fig. 22 As can be observed, in the even case 



(upper panel in Fig. 22 ) the behavior of the conductance 
is similar to a conventional NS junction with an effective 
transmission set by the second-order cotunneling ampli- 
tude. In the odd case (lower panel) the double peak 
structure within the subgap region evolves into a single 
zero bias peak as the cotunneling amplitude increases. 

To conclude this section it appears that our present 
knowledge of the non-equilibrium N-QD-S system is still 
limited and further research would be desirable, partic- 
ularly to understand the behavior of Kondo features at 
finite applied voltages and the crossover between the dif- 
ferent parameter regimes so far analyzed. In this respect 
we refer the interested reader to a recent work (Yamada 
et al , 2011| that has been published after submitting this 



review. 



C. Experimental results 

Unlike the case of the S-QD-S system, only a few works 
have addressed the issue of the transport properties of 
N-QD-S systems experimentally. This is probably due 
to the technical difficulties associated to the fabrication 



21 




FIG. 22 Non-linear conductance for t he N-QD-S system in 
the effective cotunnehng model of Ref. (Koerting et al. 2010) 
for even (upper panel) and odd (lower panel) cases. The pa- 
rameters and control the normal transmission through 
the dot in the even and odd cases respectively. The inset 
shows the differential conductance normalized to the one in 
the normal case. Reprinted figure with permission from V 
Koerting al. Physical Review B 82, 245108, 2010 ( |Koert- 
ing et al. 2010). Copyright (2010) by the American Physical 
Society. 



of such hybrid systems. The first experimental realiza- 
tion of this configuration was presented by Graber et al. 



( Graber et al. 2004 ) using a multiwall CNT as a QD cou- 



pled to Au (normal) and Al/Au (superconducting) leads 
at each side. They first analyzed the normal case by ap- 
plying a magnetic field of 25mT, clearly observing Kondo 
features in the linear conductance, as shown in the up- 



per panel of Fig. |23) At the lowest temperature of 90mK 
the normal conductance reached values ~ 1.5e^//i, indi- 
cating good and rather symmetric coupling to the leads. 
When one of the leads become superconducting it was 
observed that the temperature dependence characteris- 
tic of the Kondo regime was very much suppressed, as 
shown in the lower panel of Fig. 23 This behavior is 



in qua litative agreement with the theoretical re sults of 
Kefs. dCuevas et al] [2001[ [Tanaka et al] |2007aD for the 
case of a nearly symmetrically coupled dot. 

A different experimental realization of the N-QD-S sys- 



FIG. 23 Comparison between the linear conductance in the 
NQDN and SQDN systems as a function of the gate voltage 
for dif ferent tempera t ures i n the experimental realization of 
Ref. ( [Graber et al. 2004). Reprinted figure with permis- 
sion from M.R Grab er et al, Nanotechnology 15, S479, 2004 
(Graber 6^0^120041). Copyright (2004) by lOP Publishing 



Ltd. 



tem was presented in Refs. fDeacon et al] |2010a|b[ ). 
These authors used a self- assembled InAs QD with di- 
ameters of the order of ~ lOOnm coupled to a Ti/Au (N 
lead) and a Ti/Al (S lead). In the first of these works 
the authors focused on devices with large coupling asym- 
metry F^ ^ Fat in which the Kondo effect is suppressed 
by the strong proximity effect. In this limiting situation 
the normal lead is basically providing a means to probe 
spectroscopically the Andreev spectra of the QD-S sys- 
tem. The experiment provided evidence of the transition 
between the singlet and the doublet ground state for the 
QD-S system when the number of electrons changed from 
even to odd. As shown in Fig. |24]this was reflected in the 
behavior of the Andreev states within the gap exhibiting 
a crossing point together with a large drop in the con- 
ductance. These experimental results were in good qual- 
itative agreement with NRG calculations for the QD-S 
system. 



In a subsequent work by this group (Deacon et al. 



2010a) the same experimental realization but with vary- 



ing coupling asymmetry was analyzed. The main results 
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FIG. 24 (a) Color scale map of the differential conductance 
as a function of the source-drain and gate voltages in the ex- 
perimental realization of the N-QD-S system of Ref. (Deacon 



et al. 2010b). (b) plot of the peak conductance following the 
subgap resonances indicated by the full lines in (a). Panels 
(c) and (d) show the Andreev levels position and weight re- 
spectively obtained using the NRG method for U/A = 20 
and r^/A = 0.7. Reprinted figure with permission from R.S. 
Deacon et al. Physic al Review Letters 104, 076805, 2010 
( [Deacon et al. 2010b). Copyright (2010) by the American 
Physical Society. 



of this work are shown in Fig. [25] corresponding to asym- 
metries F^/Fat = 0.045 (upper panel) and 8.0 (two lower 
panels). For the first case with sufficiently large Tn one 
would expect the formation of a Kondo resonance due to 
the good coupling of the QD with the normal electrode. 
However, the conductance which is mediated by Andreev 
processes is suppressed inside the gap due to the very 
small coupling to the superconductor. The cases with 
asymmetries of the order of 8.0 were not in the extreme 
situation of the previous work (with F^/Fat ~ 50) and 
thus did exhibit Kondo features as can be observed in 



the lower panels of Figs. 25 



V. VOLTAGE BIASED S-QD-S SYSTEMS 

Including a finite bias voltage between the supercon- 
ducting electrodes in the S-QD-S system poses an addi- 
tional difficulty in the theory due to the intrinsic time- 
dependence of the ac Josephson effect. Even in the non- 
interacting case the inclusion of MAR processes up to 
infinite order constitutes a challenging problem which in 
general requires a numerical analysis. 

Although the mechanism of MAR processes for ex- 
plaining the subgap structure in superconducting junc- 




-0.4 -0.2 0.0 0.2 0.4 




V3, (mV) 



FIG. 25 Differential conductance at the center of odd occupa- 
tion regions for three different samples o f the experimental re- 
alization of the N-QD-S system of Ref. (Deacon et al. 2010a) 



as a function of the applied magnetic field. The three cases 
correspond to different values of the F^/Fat parameter: 0.045 
(upper panel) and 8.0 (middle and lower panels). Reprinted 
figure with permission from R.S. Deacon et al, Physical Re- 
view B 81, 121308, 2010 (Deacon et a/.| |2010at . Copyright 
(2010) by the American Physical Society. 
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FIG. 26 dc I-V characteristic for a superconducting single 
channel contact for different values of the normal transmis- 



responding to a double Fo urier transformation (Arnold[ 



1987) of G(t, t') defined as ( Martm-Rodero et al 



1999a) 



GnmM = j dt j dt'e-'^^'''-'^''^e'^^'-''^G{t,t') (36) 

The Fourier components Gnm obey an algebraic Dyson 
equation in the discrete space defined by the harmonic 
indices which can be solved using a standard recursive 
algorithm. A compact expression of these equations for 
the dot case, given in Ref. ^ Dell' Anna et al. |2QQ8 ), is 



(^Oo)nm = i^n " ^QCTz) Tz - ^ r^Cr^T^^nm 

(37) 

where Un = -\- nV, while cr^ and correspond to Pauli 
matrices in the Nambu and Keldysh space respectively 
and 



tions was introduced in the early '80s (Klapwijk et al 



1982 Octavio et a/., 1983), it was not until the mid '90s 



that a full quantitative theory of the I-V characteristics in 
a superconducting contact of arbitrary transparency was 



developed using either a scattering approach ( Averin and 

1995[) or a non-equ ilibrium 



Bardas, 1995 Bratus et al 



1996|). This x(c^) 



Green function approach (Cuevas et al 
theoretical progress allowed a very accurate desc ription 
of experimental results for atomic size contacts CScheer' 
et al] 1998a|b ). For the discussion in the present sec- 
tion it could be useful to remind the main features of 



the I-V characteristics of a one channel contact. Fig. 26 



shows the evolution of the dc current as a function of the 
contact transmission obtained using the theory of Ref. 
( Cuevas et al.\ 1996| ). As can be observed, at sufficiently 



low transmission the current exhibits a subgap structure 
with jumps aiV = 2A/n, corresponding to the threshold 
voltage for an n-order MAR process. As the transmission 
is increased the subgap structure is progressively smeared 
out and eventually at r = 1 the behavior of the I-V curve 
is almost linear except in the limit V ^ where it satu- 
rates to a finite value 2eA/h ( [Averin and Bar das 1995 



Cuevas et al. 1996) 



The case of a non-interacting resonant level coupled 
to superconducting electrodes was first analyzed in Ref. 



(Yeyati et al. 1997). We discuss briefly here the Green 



function formalism for this case which has the advantage 
of allowing to include the effect of interactions in a sec- 
ond step. The main point in this formalism is to realize 
that even when the Green functions G(t, depends on 
the two time arguments, in the case of a constant volt- 
age bias the dependence on the mean time {t^t')/2 can 
only correspond to the harmonics of the fundamental fre- 
quency eV/n dCuevas aH|1996D . This allows to express 
all quantities in terms of the components Gnm{^) cor- 



SnmX{uJn T V/2) (5n,mTl^K T V/2) \ 

Y{uj„±V/2) Sn^X{oj„±V/2) J 

(38) 

where the matrix X{uj) = —cuY/A in Keldysh space are 
given by 



2nF{uj) 
2nF{—( 



2nF{uj) 
2nF{oj) — 



\uj\ > A 
\uj\ < A 
(39) 

The presence of a discrete resonant level between the 
superconducting leads can strongly modify the I-V char- 
acteristics with respect to the quantum point contact 



case. This is illustrated in Fig. [27| which corresponds 
to a resonant level located at zero energy with decreas- 
ing tunneling rates to the leads. As the figure shows, in 
the limit of large F the I-V curves tend to that of a per- 
fect transmitting contact. In the opposite limit F <C A 
there appears a pronounced subgap structure. In con- 
trast to the contact case, the current jumps associated to 
the threshold of MAR processes appear only for the con- 
dition V = 2A/n with n being an odd integer, while the 
features at 2eA/n with even n are suppressed. This can 
be understood qualitatively from the schematic pictures 
of Fig. [28] They represent the n = 2 and n = 3 MAR 
processes with arrows indicating propagation of electrons 
(full lines) or holes (broken lines). In the n = 2 case the 
MAR "trajectory" in energy space does not cross the res- 
onant level while in the n = 3 case the resonant condition 
is fulfilled. 

The subgap features are quite sensitive to the level 
position. Fig. 29 illustrates the evolution of the I-V 



characteristics as a function of the level position eo for 
the case F^ = F^^ = A. As can be observed, when the 
level is far from the gap region the behavior of a weakly 
transmitting contact is recovered, while in the case where 
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FIG. 29 dc I-V characteristic for the non-interacting S-QD-S 
system with = = A and different dot level positions: 
5 (a), 2 (b), 1 (c), 0.5 (d) and (e) in units of A. Reprinted 
figure with permission from A. Martm-Rode ro et aL, Super- 
lattices and Microstructures 25, 925, 1999 (Martm-Rodero 



et al. 1999b). Copyright (1999) by Elsevier. 



FIG. 27 dc I-V characteristic for the no n- interacting S-QD-S 
system with eo = and different values of F = F^ = Tr. 
Reprinted figure with permission from A . Levy Yeyati et a/.. 
Physical Review B 55, 6137, 1997 (Y eyati et al.\\l997) . Copy- 
right (1997) by the American Physical Society. 





si^pressed 



FIG. 28 Schematic representation of MAR processes of order 
n = 2 and n = 3. The central horizontal line represents the 
level position, assumed to be located at the Fermi level. 



normal case. The lower panel shows the intensity map in 
the region eV < 2A for F = 0.05A. This illustrates the 
onset of higher order MAR processes, which also appear 
to be limited into wedge-like regions bounded by the con- 
dition eo = ±(A — neV/2) with odd n. The schematic 



figure 31 indicates the different resonant regions for the 
single, double and triple quasi-particle currents. 

In subsequent works further analysis of the non- 
interacting S-QD-S case out of equil ibrium was presented 
( ITonckheere a/.||2QQ9a]|Sun et a/.[[2QQ2J . While in Ref. 



Sun et aL\ 2002 ) the Hamiltonian approach was used to 



analyze the out of equilibrium dot spectral density and 
the ac comp onents of the current, in Ref. (Jonckheere] 
et al 2009a ) the effect of dephasing simulated by a third 



normal reservoir coupled to the dot has been studied. 
This work will be further commented in Sect. I VII 



A. Effect of Coulomb interactions 



the level approaches the gap, the subgap features become 
more pronounced and correspond to resonant conditions 
which depend both on A and eo- In this complex sit- 
uation a more clear picture of the overall behavior was 



presented in Ref. (Johansson et al. 1999). Figs. 30 show 



the intensity plot of the current in the cq — V plane. The 
upper panel illustrates the behavior of the current for 
eV > A for F = 0.2A, showing clearly the onset of single 
quasiparticle current for eV > 2 A at eo = 0. For eo 7^ 
this current is only significant in a wedge-like zone lim- 
ited by eo = ±(V — 2A)/2. It can also be noticed in 
addition the presence of resonant peaks at V/2 = ±eo 
which are reminiscent of the resonant condition for the 



The inclusion of intradot interactions in the out of 
equilibrium S-QD-S system introduces an additional dif- 
ficulty in an already challenging theoretical problem, as 
shown in the previous section. So far the attempts have 
been restricted to some limiting cases which have been 
treated using approximate methods. One of these spe- 
cial limiting situations which was first analyzed was the 
case of a quantum dot in the strong Coulomb blockade 



regime (Whan and Orlando, 1996 Yeyati et al] 1997). 



These works were motivated by the experimental results 



of Ref. (Ralph et al. 1995) for transport through small 



metallic nanoparticles coupled to superconducting leads. 
In this strong blockade regime multiple quasiparticle pro- 
cesses are suppressed and the current is basically due to 
single quasiparticle tunneling. 
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FIG. 31 Lines indicating the boundaries for the resonant re- 
gions of the single, double and triple quasiparticle currents 
in the eo-eV plane for the non-interacting S-QD-S system. 
Reprinted figure with permission fr om G. Johansson et aL , 
Physical Review B 60, 1382, 1999 ([Johansson et al.\\1999\. 



FIG. 30 Color scale map of the the dc current in the eo- 
eV plane for the non-interacting S-QD-S system. The upper 
panel corresponds to eV > A and F = 0.2 A while the lower 
panel corresponds to eV < 2A and F = 0.05A. Reprinted 
figure with permission fr om G. Johansson et aL, Physical Re- 
view B 60, 1382, 1999 ( [Johansson et al.\ 11999] ). Copyright 
(1999) by the American Physical Society. 



Copyright (1999) by the American Physical Society. 



pie sequential tunneling picture, which would predict 
I{V) - rL{e)rR{e)/{rL{e) + r^(6)), exhibiting an in- 
trinsic broadening of the BCS-like feature, in agreement 
with the experimental observation (Ralph et al.\ 1995). 



A similar expression was obtained in Ref. ( [Kang 1998) 



In Ref. (Whan and Orlando, 1996) the current was jj^ j^^f (jAvishai et al. 



calculated in this regime by means of a master equation 
approach assuming a sequential tunneling regime. The 
single-particle tunneling rates were calculated using the 
Fermi go lden-rule. A s lightly different method was used 
in Ref. ( Yeyati et aL\ ^97) where resonant tunneling 
through an effective one-electron level describing the dot 
in the limit U ^ oo was considered. The corresponding 
expression for the tunneling current was given by 



using the equation of motion approach in the atomic 
limit which produces a correction factor in the current, 
~ "^^{l— < TiQcr >) due to the strong Coulomb interac- 
tion. 

In order to analyze the interplay of MAR and Kondo 
correlati ons it is necess ary to rely on other approaches. 

2003) the current- volt age of this 



system was obtained using the slave boson mean field 
approximation already discussed in Sect 



III 



In the 
([AvishaiJ 

et al] [2003|) the problem becomes equivalent to an effec- 



infinite-U version of the method used in Ref. 



tive non-interacting model with renor malize d parameters 

The authors 



III.B 



4e 
~h 



I 



duj- 



(c^-e)2 + [r^(c^)+r^(a;) 



X [npiuj - eV/2) - npiou + eV/2)] , (40) 
where e denotes the effective level and 
TlM^) = (r/2)Re [|c^ ± F/2|/V(^±F/2)2- A2 . 

The corresponding result for different values of F are 
shown in Fig. [32] This result differs from the sim- 



F and eo, as indicated in subsection 
analyze the evolution of the I-V curves and the shot-noise 
as a function of the parameter T^/A. As is physically 
expected, the behavior is similar to that of a perfectly 
transparent contact when T^/A 1 developing a clear 

subgap structure in the opposite limit. 

A similar approach was used in Ref. (Yeyati et al 



2003) in order to analyze the low bias transport proper- 



ties of a S-QD-S in the Kondo regime Tx ^ A. It was 
shown that these properties can be understood in terms 
of the dynamics of the subgap Andreev states. In this 
limit the ABSs satisfy the equation corresponding to the 
non-interacting case, ie. Eq. (12), with renormalized pa- 
rameters Tl^r and eo instead of Tl^r and eo- The low 
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FIG. 32 dc current-voltage characteristic for a S-QD-S system 
in the strong Coulomb blockade regime for different values of 



= 5 X IQ-'^A (full line), IQ-'^A (dashed line) and 2 x 
10~^A (dotted line), with e = 5A. Reprinted figure with 
permission fro m A. Levy Yeyati et aL, Physical Review B 
55, 6137, 1997 ([Yeyati et a/.1[T997t. Copyright (1997) by the 



American Physical Society. 



bias quasiparticle current through the system arises from 
transitions between the continuum occupied and empty 
states below and above the superconducting gap which 
occur by means of Landau- Zener like processes involving 
the ABSs. This is illustrated in Fig. [33| first for the case 
of a quantum point contact (left panel) and then for the 
S-QD-S case (right panel). In the last case it is necessary 
to have a transition between the ABSs and the contin- 
uum in addition to the Landau-Zener transition between 
the lower and upper ABSs. The results for the low bias 
dc current which are obtain from this analysis, shown 



in Fig. 34, are in good agreement with the results of a 
full numerical calculation including MAR processes up 
to infinite order. The analysis based on the dynamics 
of the ABSs was also used in Ref. (Vecino et at., 2004) 



for a comparison of experimental results in the low bias 
regime. In contrast to Ref. (Yeyati et al, 2003), in this 
work a phenomenological damping rate t] was introduced 
in order to fit the experimental data obtained for multi- 
wall CNTs connected to Au/Al leads. 

A step beyo nd the infinite- SB MF approach was 
taken in Ref. (Eichler et al, 2007) where the finite-U 



SBMF method was used to determine the I-V charac- 
teristics of a S-QD-S system. This approach allowed to 
describe the observed differences in the subgap structure 
between situations with even and odd number of elec- 
trons in the dot for a SWCNT QD c oupled to Al/Ti 
electrodes. In Ref. (Eichler et al., 2007) the effective pa- 
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FIG. 33 Schematic representation of the ABS dynamics ex- 
plaining the low bias dc current in the case of a single channel 
contact (panel a) and for the S-QD-S system in the Kondo 
regime (panel b). Reprinted figure with permission from A. 
Levy Yeyati et aL, Physical Review Letters 91, 266802, 2003 
( Yeyati et aZ"||2003| ). Copyright (2003) by the American Phys- 
ical Society. 
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FIG. 34 Comparison of the low bias dc current in the S- 
QD-S system in the Kondo regime obtained from the ABS 
dynamics (dashed lines) and by the full numerical calculation 
(full line) for different values of the effective coupling to the 
leads Tl — — 1, 2, 3, 4, 5 and 10 in units of A, from bottom 
to top. Reprinted figure with permission from A. Levy Yeyati 



et al Physical Review Letters 91, 266802, 2003 (j Yeyati et al. 
2003). Copyright (2003) by the American Physical Society. 



rameters of the finite-U SBMF approach were obtained 



for the leads in the normal state and assumed to remain 
unmodified in the superconducting case. The results of 
this work for the subgap structure are shown in Fig. 
35 where the comparison with the experimental data is 
given. The figure corresponds to an odd valley exhibit- 
ing clear Kondo features in the normal state. The more 
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FIG. 35 Experimental and theoretical results for the differ- 
ential conductance of the S-QD-S system in an odd valley. 
Panel (a) corresponds to the normal state (experimental and 
theoretical results), while (b) and (c) show the correspond- 
ing results in the superconducting state. Panel (d) shows the 
conductance vs voltage bias along the line indicated in red 
in panels (b) and (c). The theoretical curves corresponds to 
Fl = a, Fr = 0.03A and U = lOA. Reprinted figure with 
permission from A. Eichler et aL, Physical Review Letters 99, 

Copyright (2007) by the 



2007) 



126602, 2007 (Eic hler et al. 
American Physical Society. 



intriguing feature of the differential conductance in the 
superconducting case was the presence of a pronounced 
structure for eV ~ A, which cannot be explained by 
the non-interacting theory. This feature was attributed 
to the large asymmetry Fl/Fr ~ 30 which produces a 
Kondo resonance pinned at the chemical potential of the 
left lead. As is schematically depicted in Fig. |35|e) this 
resonance would produce and enhancement of the current 
for eV ~ A. 

The effect of interactions in t he subgap structure has 
also been considered in Ref. (Dell' Anna et al. 2008) 



by means of a perturbative approach in which the dot 
self-energy was calculated up to second order in U. The 
method thus includes the diagrams already discussed in 
Sect. Ill but extended to the non-equilibrium situa- 



tion. To avoid heavily time-consuming computation of 
the multiple frequency integrals, the authors calculate 
the diagrams in time representation and then Fourier 
transform the final result. They consider the weak in- 
teraction regime U/F < 1 thus avoiding the regime of 
TT-j unction behavior. The most remarkable result is the 
observation of an enhancement of the current due to the 
interactions, which is more pronounced for voltages ap- 
proaching the odd MAR onset conditions 2A/(2r + 1). 
This is illustrated in Fig. [36] for the case eo = where 
the difference I{U) — /(O) as a function of 2A/V is rep- 
resented. The enhancement is observed both in the self- 
consistent first order approximation and when including 
the second-order diagrams. This current enhancement is 
reminiscent yet different from the "antiblockade" behav- 
ior due to dynamical Coulomb blockade effects on MAR 




FIG. 36 Variation of the dc current in the symmetric S-QD- 
S syste m with respect to the n on-interacting case obtained 
in Ref. ( Dell' Anna et al.\ 2008) using the second-order self- 
energy approach for U/F — F/A = 0.5. The dashed curve 
gives the self-consistent first order result. Reprinted figure 
with permission from L. Dell'Anna et al, Physical Review B 
77, 104525, 2008 (DelPAnna et al.\^OOS\. Copyright (2008) 



by the American Physical Society. 



transport, as discussed in Ref. (Yeyati et al , 2005). It is 



worth mentioning that MAR transport through a reso- 
nant level co upled to a localized pho non mode was stud- 
ied in Ref. ( Zazunov et al. , 2006 ) using second order 



perturbation theory in the electron-phonon coupling. 



B. Summary of Experimental results 



et 



The alrea dy commented work by Ralph et al. (Ralph 
T995| can be considered one of the first realizations 



of a S-QD-S system in which the current-voltage charac- 
teristic was measured. This case corresponded, however, 
to the strong blockade regime in which the subgap struc- 
ture is absent. It was not until 2002 that experiments on 
CNTs coupled to Al leads (Buitelaar et al.^ ^2003, 2002] ) 
allowed a clear observa tion of the subgap features. Ref. 
(Buitelaar et al. 2002 ) mainly focused in the linear con- 



ductance which can exhibit either an enhancement or a 
suppression with respect to the normal case depending 
on the ratio TWA. The results of this work are summa- 
rized in Fig. 37 On the other hand in Ref. (Buitelaar 
(eTa/.l |2QQ3) the authors analyzed the MAR induced sub- 
gap structure for the same type of systems in more detail. 



As shown in Fig. 38 clear peaks in the differential con- 
ductance are observed at the positions eV ~ 2 A, A and 
A/2. The intensity of these peaks evolves as a function of 
the dot level position (controlled by the gate voltage Vg). 
Contrary to the theoretical expectations the peak at A 
is still visible at resonance whereas the expected feature 
from the non-interac ting model at 2 A/3 is not observed. 
The authors of Ref. ( [Buitelaar et al. 2003 ) suggest that 
the discrepancy can be attributed to the effect of inter- 
actions not included in their theoretical analysis. These 
results could be analyzed in the light of the already com- 
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FIG. 37 Linear conductance in the Kondo regime for the ex- 
perimental reahzation of the S-QD-S system of Ref. ( Buite- 
laar et al. 2002| normahzed to its value in normal state as 
a function of Tk/^- Reprinted figure with permission from 
M.R. Buitelaar et a/. , Phys ical Review Letters 89, 256801, 
2002 (Buitela ar et a/.[|2002t . Copyright (2002) by the Amer- 
ican Physical Society. 



mented arguments of Ref. (Eichler et al. 2007) which 
attributed the pronounced A peak to the combined ef- 
fect of coupling asymmetry and Kondo effect. 

The competition between Kondo ef fect and supercon- 
ductivity was also observed in Ref. (Buizert et al, 2007) 



in which self-assembled InAs quantum dots contacted 
with Al leads were analyzed. A magnetic field was used 
to control the size of the superconducting gap parame- 
ter and the linear conductance as a function of A/ksTK 
was measured. The results exhibited a rather univer- 
sal behavior as a function of this parameter. However, in 
contrast to Ref. (|Buitelaar et al , 20021) the ratio Gs/Gn 
did not exceed unity for A/Zc^T^ < 1 while decreasing 
as expected for A/kBTx > 1. Although the authors at- 
tributed this difference to a stronger Coulomb repulsion 
which would in their case heavily damp the MAR pro- 
cesses, one would expect that this effect would be already 
included when scaling A in units of Tk- The absence 
of conductance enhancement for large Tk could be also 
pointing out to an ingredient in this system not included 
in the simplest Anderson model like spin-orbit interac- 
tions. 



VI. BEYOND THE SINGLE LEVEL MODEL: MULTIDOT, 
MULTILEVEL AND MULTITERMINAL SYSTEMS 

In recent years there has been an increasing interest in 
more complex situations which cannot be described by 
the simplest single level Anderson model. These include 
situations where transport occurs through more than a 
single dot or where several quantum channels in a sin- 
gle dot are involved. In addition, there is also great 
interest in analyzing the transport properties of hybrid 
quantum dot systems coupled to several superconducting 
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FIG. 38 Results for the differential conductance in the exper- 
imental realization of the S-QD-S system of Ref. ( Buitelaar] 
let g/., 2003). Panel (a) color map for this quantity in the bias 
voltage-gate voltage plane. The dashed lines indicate the evo- 
lution of the MAR resonances. Panels (b) and (c) show plot of 
the conductance as a function of bias for the lines indicated 
by the arrows in panel (a). Reprinted figure with permis- 
sion from M.R. Buitelaar et a/., Physical Review Letters 91, 
057005, 2003 (B uitelaar et al] |2002t . Copyright (2003) by 
the American Physical Society. 



and/or normal electrodes in a multiterminal configura- 
tion. These configurations could allow to explore non- 
local electronic transport, in particular the possibility of 
creating entangled electron pairs by means of crossed or 
non-local Andreev processes. To describe these devel- 
opments we organize this section as follows: in the first 
subsection we discuss the case of Josephson transport 
through double and multiple dot systems, in the second 
one we consider this effect for a multilevel dot, and finally 
we consider multiple dot systems including both normal 
and superconducting electrodes as well as in a multiter- 
minal configuration. 



A. Josephson effect through multidot systems 

Transport through double quantum dots connected ei- 
ther in parallel or in series to superconducting electrodes 
has been extensively analyzed in the literature. Most 
of the theoretical works describe this situation by us- 
ing a single-level Anderson model to represent each dot 
and introducing extra terms describing the coupling to 
the leads. Choi et al. considered in Ref. (Choi et al. 



2000) the case of two dots connected in parallel to su- 
perconducting leads as depicted in Fig. [39] By ana- 
lyzing the problem to the fourth order in the tunneling 
to the leads they derived an effective Hamiltonian cou- 
pling the localized spins in both dots. In the regime 
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FIG. 39 Schematic representation of a double dot system 
coupled in parallel to two superconducting electrodes (upper 
panel) considered in Ref . ( Choi et al. 2000 ) . The lower panel 



describe the leads spectral density and the parameters of the 
double dot Anderson model used in this reference. Reprinted 
figure with permission from M.S. Choi et aL, Physical Review 
B 62, 13569, 2000 ( |Choi et a/.[[2000| . Copyright (2000) by 
the American Physical Society. 



< —€q « A « U this Hamiltonian adopts the form 



Hef / J (1 + COS (f) 



Sa-Sfo — - 



(41) 



where J is an exchange coupling between the localized 
spins in dots a and 6; and cp is defined as 



if = (1)1 



dlb) .A, 



(42) 



where the last term corresponds to the phase accumu- 
lated on each path of the loop due to the magnetic 
field. These results indicate that the Josephson current 
through such a device would be sensitive to the total 
spin of the double dot. In order to probe the spin state 
of the double dot system the authors propose to incor- 
porate it into a SQUID geometry in which an additional 
tunnel junction is included in one of its arms. Further 
elaboration on similar ideas were presented in Ref. (Hur 
\et al.\ [2006] ) for a triple dot between superconducting 
leads. This system is shown to behave under certain 
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FIG. 40 [/ — e phase diagrams, dot population and spin cor- 
relation functions for the serie s double dot model considered 
in Ref. (Bergeret et al. 2006| ) obtained using the zero band- 
width approximation to describe the superconducting leads. 
The parameters are ti2 = 10 A, ti, = = 2 A (a) and 
th — tn — 2. 5 A (b). The plots in panel (c) are taken along the 
dotted line corresponding to [/ = 800 A shown in panel (a). 
Reprinted figure with permission from F.S. Bergeret et al , 
Physical Review B 74, 132505, 2006 ( |Bergeret et a/.[|2006| ). 
Copyright (2006) by the American Physical Society. 



conditions as a mesoscopic pendulum where the singlets 
injected through a pair of dots oscillate between two dif- 
ferent configurations like in the resonating valence bond 
model. 

The case of two QDs in series conne cted to supercon - 
ducting leads was first analyzed in Ref. ( Zhu et al] 2002 ). 



The authors start by diagonalizing exactly the isolated 
double dot Hamiltonian including a term describing the 
interdot repulsion V and then introducing the external 
leads by a Dyson-type equation. This approach should 
be valid in the limit of vanishing coupling to the leads 



and is equivalent to the one commented in Sect. Ill for 
the single dot case. For the case U > V > where t is 
the interdot hopping parameter, and for ksT F they 
find that the system exhibits a type current-phase re- 
lation except when a finite Zeeman splitting is included. 
The absence of a 7r-j unction behavior for any dot filling is 
probably due to the extremely small value of F compared 
to ksT. This will be further discussed below. 

A similar situation was later considered in Ref. (Berg- 



eret et al] 2006} , where a double dot model was analyzed 
using both small cluster numerical diagonalizations (dis- 
cussed before in Subsect. III.D| for the single dot case) 
together with the finite-U SBMF technique. In addition 
to the transition to the 7r-phase this work aimed to inves- 
tigate the interplay between different possible magnetic 
correlations including Kondo and anti-ferromagnetic cou- 
pling between the localized spins within each dot. Part of 



the results are shown in Fig. 40 where the phase-diagram 
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FIG. 41 
in Ref. (Bergeret et al 



Current-phase relation for the model considered 
06) obtained using the finite- 
U SBMF method with parameters U = 800A, ti2 = lOA, 
Tl =Tr = 2.25A (upper panel) and Tl = = 4A (lower 
panel). Reprinted figure with permission fro m F.S. Bergeret 
et al, Physical Review B 74, 132505, 2006 ([Bergeret efal\ 



2006). Copyright (2006) by the American Physical Society. 



in the U vs — e plane is given for the case ti2 ^ A, ti2 
being the interdot hopping parameter. In the range of 



parameters of Fig. 40 the system exhibits a 7r-phase re- 
gion associated to the transition between the empty and 
singly occupied double dot. When increasing the dot 
population to the level of one electron per dot antifer- 
romagnetic correlations dominate and the 7r-phase is no 
longer stable. This is further illustrated in panel (c) of 
Fig. 40 where different spin-spin correlation functions 
are shown both for the normal and the superconducting 
case. The results obtained using the small cluster diago- 
nalizations were confirmed by the finite-U SBMF calcula- 
tions. In particular these last calculations also show the 
gradual disappearance of the full 7r-phase when increas- 



ing the coupling to the leads, as illustrated in Fig. 41 



This work also provided a possible scenario for explaining 



the experimental results of Ref. ( jKasumov et al\ [2005) 



for fullerene dimers containing Gd magnetic atoms and 
suggested the possibility to control the magnetic configu- 
ration of these atoms by means of the Josephson current. 

The SBMF approach in the infinite-U version was ap- 
plied to series and p arallel double quantum dots in Ref. 
(Lopez et al. 2007). As in the case of a single QD 



this method cannot account for the appearance of a tt- 
junction phase. The results are nevertheless relevant for 
the regime Tk ^ A where Kondo correlations dominate 
over pairing. While in the parallel case it is found that 
the Josephson critical current, /c, decreases monotoni- 
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FIG. 42 Josephson critical current as a fu nction of t/F for the 
series double dot model analyzed in Ref. (Lopez et al. 2007) 
using the infinite-U SBMF method for A/Tr 0.1,0.25 and 
0.5 (from bottom to top). The upper inset shows the compar- 
ison of the numerical result for A/Tk =0.1 (full line) with 
the prediction from a non-interacting model with an effective 
transmission (dashed line). The lower inset shows the criti- 
cal phase at which the maximum current occurs. Reprinted 
figure with permissio n from R. Lopez e t a/., Physical Review 
B 75, 045132, 2007 ( Lopez et a/.][2007| . Copyright (2007) by 
the American Physical Society. 



cally with the interdot hopping parameter in the se- 
ries case a non-monotonous behavior is found. This is 
illustrated in Fig. 42 where Ic exhibits a maximum at 
t/r ~ 1. The authors interpret the change in behav- 
ior of Ic as a transition from a regime characterized by 
two independent Kondo singlets involving each dot and 
the corresponding lead to the formation of bonding and 
antibonding Kondo resonances. 

Quite recently, the series double dot system coupled 
to SC leads has be en analyzed using the NRG method 
(Zitko et al. 2010). The authors considered the regime 
U ^ oo and —e:^r, which would correspond to the 
deep Kondo regime in a normal single QD. In this range 
of parameters they find a rich phase diagram as a func- 
tion of the interdot coupling t and the ratio T^/A. Some 
of their results are illustrated in Fig. [43] showing the re- 
gions corresponding to 0, tt and tt' phases in the t/T vs 
plane for different values of T/^/A. One can notice the 
abrupt transition between and tt phases for t ^ lOF 
which can be associated to a change in the DQD popu- 
lation from an even to an odd number of electrons. An 
additional remarkable feature is the appearance of a tt' 
"island" close to ^ = tt and for t ~ F in the intermediate 
coupling regime Tk ~ A. 

The evolution of the Josephson effect as a function of 
the number of dots connected in series was studied in 
Ref. 



(Bergeret et al. 



2007). The model considered in 



this work with all dot levels fixed at a same vale e, with 
the same local Coulomb repulsion U and with dots con- 
nected by a fixed hopping parameter is equivalent to a 
finite Hubbard chain connected to two superconducting 
leads. The ground state properties of this model were 



31 



r=0.014D 



r=0.02 D 



c) 

1 Uo- 



10" 

10"" ^ 
10' 
10^ 
10"' 
10"^ 



a) 



I I I 



I I I 



I I I 



b) 



Q GQJD O- G G 



I I I I I . I 



Ti: 



■o o- 





I I I I I I I 



-O Q G- 





I I I I I . I 



Strong 
coupling 



A/T^ =3.6 

Intermediate 
coupling 



A/Tj^ =10 

Weak 
coupling 



0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 1 

(t)/7U 

FIG. 43 Phase boundaries betwe en the and tt sta tes in the 
series double dot model of Ref. (IZitko et al.\ 12010') obtained 



using the NRG method for F = 0.014^ and F = 0.02^, where 
D is the leads band-width. Reprinted figure with permission 
from R. Zitko et aL, Physical Review Letters 105, 116803, 
2010 dZitko et a/.]|2010| ). Copyright (2010) by the American 
Physical Society. 



obtained using the zero band-width limit description of 
the leads discussed in Subsect. IILD and employing the 
Lanczos algorithm. Fig. 44 illustrates the evolution of 



the phase diagram as the number of dots in the chain, 
is increased. One can clearly distinguish the case of even 
and odd TV. In the last case the diagram is similar to the 
single dot case with a central 7r-phase region correspond- 
ing to the half- filled case. On the contrary, for even N 
the TT-phase is absent around half-filling due to the dom- 
inance of antiferromagnetic correlations between spins in 
neighboring dots. One can also notice the appearance 
of additional narrower regions of 7r-phase character cor- 
responding to fillings with odd number of electrons in 
the dots region. The authors also analyzed the current- 
phase relation as a function of N for the half- filled case. 
Fig. 



45 shows that the critical current scales as e" 



-aN 



with a different sign depending on the parity of N. This 
behavior is consistent with the prediction of field theoret- 
ical calculations for a ID Luttinger liquid with repulsive 
interactions where the fixed point corresponds to the ab- 
sence of Jos ephson coupling in the limit of an infinite 
long chain Affleck et a/.[|2000 ). 
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FIG. 44 e — U phase diagrams for the quantum d ot ar- 
ray coupled to superconductors considered in Ref. ( |Berg- 



eret et al. 2007) with increasing number of dots N = 2,3,4 



and 5. The results were obtained using the zero band- width 
approximation described in Subsect. |III.D| with parameters 
tL = tR = t = A. Reprinted figure with permission fro m F.S. 
Bergeret et al, Physical Review B 76, 174510, 2007 (iBerg- 



eret et a/.[|2007t . Copyright (2007) by the American Physical 
Society. 
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FIG. 45 Josephson current-phase relation for the qua ntum 
dot array model analyzed in Ref. ( [Bergeret et al.\ [200 7^ with 
U — lOA at half-filling. The curves in the positive half plane 
correspond to A/" = 2, 4 and 6 (from top to bottom) while those 
taking negative values correspond to A^ = 1,3 and 5 (from 
bottom to top). The current is scaled by and exponential 
factor exp aN with a ^ 1.8 in units of eA/h. Reprinted figure 
with permission from F.S. Bergeret et al, Physical Review B 
76, 174510, 2007 (B ergeret et a/.|[2007| . Copyright (2007) by 
the American Physical Society. 
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FIG. 46 Josephson current at — 7r/2 in the multilevel S- 
QD-S model of Ref. ( [Shimizu et al. 1998| as a function of the 
level position for the cases n = 1 (a), n = 2 (b) and n = 3 (c), 
where n denotes the level degeneracy. The dotted lines indi- 
cate the evolution of the dot total charge Ud- Reprinted figure 
with permission from Y. Shimizu et a/., Journal of the Phys- 
ical Society of Japan 67, 1525, 1998 ([Shimizu et~ar| |1998 ). 
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B. Multilevel quantum dots 

So far only the single-level Anderson model has been 
considered for describing an individual quantum dot. A 
proper description of actual physical realizations of quan- 
tum dots could require to consider a multilevel general- 
ization of this model. This has been already p ointed out 
in connection with the experiments of Ref. (van Dam 



et fl/.[|2QQ6J on In As nanowires, whose results where qual- 
itatively accounted for using a multilevel model in which 
two orbitals with different parity were involved. 

Multilevel effects in the Jose phson current throug h a 

by 



1998 



QD were first addressed in Ref. ( Shimizu et al , 
means of the Hartree-Fock approximation. The authors 
showed that when non-diagonal processes involving dif- 
ferent dot levels are relevant the system can behave as 
a TT-j unction even in the absence of a magnetic ground 
state. This is illustrated in Fig. [46] where the Joseph- 
son current at (j) = 7r/2 is plotted as a function of the 
dot levels position, eo, for the cases with 1, 2 and 3 lev- 
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FIG. 47 Phase dia gram in the J — A plane for the two level 
dot model of Ref. (Lee et al. 2010') corresponding to an ef- 



fective one channel situation with symmetric coupling to the 
leads. The notation for the different phases is indicated in 
the main text. Reprinted figure with permission from M. Lee 
et a/.. Physical Review B 81, 155114, 2010 ( |Lee et a/.[[2010 ). 
Copyright (2010) by the American Physical Society. 



els and including diagonal and non-diagonal couplings to 
the leads with the same value F. One can clearly notice 
that whereas in the lower panel (single level) the 7r-phase 
is only present for odd number of electrons {rid = 1) in 
the 2 and 3 level cases (panels b and c) the tt behavior is 
also present for even occupancy plateaus. 

Similar ideas were discussed in Ref. (jRozhkov et al. 



2001) for a multilevel situation with nearly degenerate 
levels (^e <C A) connected at the same two points to the 
leads and using the cotunneling approach. The authors 
pointed out that the 7r-j unction behavior is linked to 
two-particle processes in which one of the electrons pro- 
ceeds through an occupied state and the other through 
an empty one. 

Furthe r analysis of this multilevel case was presented 
in Ref. ( Lee et aL\ 2010[ ) using both perturbation the- 
ory (cotunneling approach) and NRG calculations. They 
considered a two-level situation coupled in parallel to 
single channel leads and including an exchange term J 
between the electron spins in each dot level. In the 
normal case this exchange term for J < would drive 
a singlet-triplet transition. Two different situations are 
distinguished: a case in which the two orbitals have the 
same "parity" (i.e. sign(tiL * tm) = sign(t2L * hn), 
where tj^ are the hopping from the j level to the lead 
a) and a case in which the parities are different. In the 
first case and for tiL/tiR = ^2L/^2i? the dot levels are 
only coupled to the symmetric combination of the two 
leads, yielding an effective one channel problem. In the 
second case the problem is equivalent to a two channel 
two impurity model with exchange coupling. 
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FIG. 48 Same as Fig. [47| for the case corresponding to an 
effective two-channel two-impurity model. Reprinted figure 
with permission from M. Lee et al, Physical Review B 81, 
155114, 2010 dLee et aL\ |20To| ). Copyright (2010) by the 
American Physical Society. 



The system exhibits a very rich phase diagram depend- 
ing on the several model parameters. For illustration we 
show in Fig. 



47| and [48| the obtained phase diagrams in 
the A, J plane for the single channel and the two chan- 
nel cases respectively for symmetric coupling to the leads. 
The different phases are denoted by a capital letter indi- 
cating the spin of the dominant ground state, which can 
be either 5, or T for spin 0, 1/2 or 1. In the case 
of mixed phases where the ground state changes with 
(j) a subindex is included indicating the character of the 
metastable state. There is finally a label which can be ei- 
ther or TT indicating the character of the current-phase 
relation in the dominant phase. The different types of 
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current-phase relations are shown as insets in Figs, 
and SHI 

Another multilevel situation which has been recently 
analyzed is the case of a fo ur-fold degenerate lev el cou- 
pled to multichannel leads (Zazunov et al., 2Q1Q| ). This 



situation has been found in experiments on high qual- 
ity CNT quantum dots with normal leads showing clear 
signatures of four-fold degeneracy both in the Coulomb 



blockade and in the Kondo regimes ( Jarillo-Herrero et al 



2005 2007). The SU(4) Kondo effect has been analyzec 



1999| ) , dChoi et al 
Ref. ( [Zazunov et al 



in the normal case by several authors like (Yeyati et al. 

and ( |Lim et al] |2QQ6D . 



2005b) and ( |Lim et al] |2QQ6D . In 
2010^ Zazunov et al. have studied 



the Josephson effect in this case by considering a general- 
ized SU(4) Anderson model with superconducting leads. 
They obtained analytical results in two opposite regimes 
corresponding to the deep Kondo limit Tk ^ A and 
the cotunneling limit. In the first case they obtained a 
current-phase relation which corresponds to the super- 
position of effective non-interacting channels with renor- 




FIG. 49 Phase diagram in the U — e plane for the SU(4) 
Anders on model with superconducting leads analyzed in Ref. 
( Zazunov et al. 2010) obtained in the ZBW approximation 
for A = lOr. White regions correspond to the (S, T) — (0, 0) 
and green to {S,T) = (1/2, 1/2). The black and blue regions 
correspond to mixed phases with (S, T) = (0, 0) at = and 
(5',T) = (1/2,1/2 or (5', T) = (1,0) - (0,1) at = tt respec- 
tively. The right panel shows the evolution of the current- 
phase relation for e/F = —5 and several values of U/T. The 
small region in red inside the dashed white line corresponds to 
a 7t' mixed phase. Reprinted figure with permission from A. 



Zazunov et al, Physical Review B 81, 012502, 2010 (Zazunov 



et al] 2010). Copyright (2010) by the American Physical So- 



ciety. 



malized transmission r = 1/2, therefore deviating from 
the SU(2) case where the unitary limit can be reached. 
In the cotunneling case with U ^ oc a O-tt transition at 



was obtained, as i n the SU(2) case (see Ref. ( Glaz 

"man and Matveev' "1989')) but with a different ratio of the 
critical currents Ic{—\eo\)/Ic{\eo\) = —1/4. The reduc- 
tion of this ratio by a factor 2 with respect to the SU(2) 
case can be readily understood by considering the num- 
ber of processes leading to Cooper pair transfer through 
the dot. 

The authors also analyzed numerically the phase di- 
agram of the model in the regime A F. A first 
insight is obtained by taking the limit A ^ oo in 
which case the relevant Hilbert space is reduced to the 
2^ dot states. Conservation of the total spin, and 
orbital pseudo-spin, T, allows to further decouple this 
Hilbert space into three different sectors: (S^T) = (0,0), 
{S,T) = (1/2,1/2) and {S,T) = (1,0) or (0,1) (these 
last two are degenerate in the SU(4) case). The main 
features of the phase diagram obtained in this limit were 
shown to be preserved in the finite A ^ F regime, which 
was analyzed by means of the zero band-width model for 
the leads. The phase diagram, illustrated in Fig. [49) 
is essentially the same as in the A ^ oo limit except 
for two properties: first, the appearance of tiny tt' type 
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mixed phase (indicated by the red regions in Fig. 49) 



and second by the change in the character of the current- 
phase relation of the {S^T) = (1/2, 1/2) with increasing 
U. The panel on the right shows that this relation evolve 
with U/A: while for A U it is typically of type, for 
U > A it becomes of tt type. 

In closing this subsection we mention the recent ap- 
pearance of a work by Lim et al. analyzing the effe ct of 
including spin-orbit interactions within this model (Lim 
et aL\ 201 1[), which can be relevant for small radius CTNs 



(Kuemmeth et a/., 2008) 



C. Multidot-multiterminal systems with normal and 
superconducting leads 

1. Josephson effect through a quantum dot in a three terminal 
configuration 

The Josephson effect through a single dot coupled to 
two superconductors and to a third normal lead has 
been analyzed by several authors. This configuration is 
schematically depicted in the upper panel of Fig. [50] An 
interesting non-equilibrium enhancement of the Joseph- 
son effect was predicted in Ref . ( Pal a et al] |2007| ) . The 
authors applied the real time diagrammatic approach 
commented in Subsect. |III.C| with the tunneling rates 
to the leads calculated to the first order in Ts^n- They 
found that a significant Josephson current can be induced 
by the voltage bias applied to the normal leads even in 
the situation F^ < ksT where the equilibrium Josephson 
current would be negligible. 

Fig. 50 shows a color map of the Josephson critical 



current as a function of both the level position and the 
chemical potential on the normal lead, jun- One can no- 
tice the presence of a white region for \jiiN\ U/2 where 
the Josephson current is negligible. Outside this region 
it becomes of the order of ^ eTs/h and exhibits a tran- 
sition from to TT behavior depending on the level posi- 
tion. The origin of this peculiar behavior can be traced 
to the enhancement of the proximity effect pairing ampli- 
tude on the dot due to the non-equilibrium population 
which increases the double population probability that 
is strongly suppressed at /iAr = due to the charging 
energy. In a sub sequent publication by the same group 
( Governale et a/.^ 2008) the authors considered the same 
effect in the limit A ^ oo which allows to account for 
the Josephson effect to all orders in F^. The results ob- 
tained are qualitatively similar and allow to identify the 
lines separating different regions with the Andreev bound 
states of the S-QD-S system. 

It should be also mentioned within this context the 
work by Jonckheere et al. fJonckheere et a/., 2009b) in 
which the effect of a third normal lead on the ac Joseph- 
son effect in a non-interacting voltage biased S-QD-S sys- 
tem was analyzed. The main idea of this work was to 
study the transition from the coherent MAR regime to 
the incoherent limit controlled by the coupling F^r to the 
normal lead. They show that while the dc Josephson 




FIG. 50 Color map of the Josephson critical current (lower 
left panel) and Andreev current (lower right panel) for the QD 
dot coupled to two superconducting leads and an additional 
voltage biased normal one in the iin /U — e/U plane. The up- 
per panel gives a schematic representation of the setup con- 
sidered in this work. Reprinted figure with permission from 
M. Pala et al, New Journal of Physics 9, 278, 2007 ( |Pala| 
eran[2007|. Copyright (2007) by lOP Publishing Ltd. 



current exhibits a monotonous decrease with increasing 
Ftv the behavior of the dc quasiparticle current and its 
first ac harmonics have a much more involved evolution. 



which is illustrated in Fig. 51 



2. Andreev transport through double quantum dots 

Andreev transport in double quantum dot systems con- 
nected to normal and superconducting leads has been 
studied so far in a few works. Tanaka et al. ( jTanaka 



et al. 2008) considered the case of a T-shape geometry 



where a central dot is coupled to both electrodes and 
a second dot is only side-coupled to the central one, as 
shown in the upper panel of Fig. [52j The authors used 
the NRG method in the A ^ oo limit where, as discussed 
in Subsect. |III.D[ the superconducting lead acts as a sim- 
ple boundary condition for Andreev reflection. They fo- 
cused in the case where interactions are neglected in the 
central QD analyzing the effect of increasing U in the 
lateral dot on the conductance through the system. The 
results of Fig. [52] show that the Andreev conductance 
gradually approaches the unitary limit as the side dot U 
increases for the case F^v = where t is the hopping be- 
tween the two dots. On the other hand, for smaller values 
of Fat the evolution of the conductance with U is non- 
monotonous exhibiting for the symmetric case first an 
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FIG. 51 dc and the two first harmonics of the current- voltage 
characteristic of a non-interacting single level S-QD-S system 
with an additional normal lead connected to the dot with 



increasing coupling Tn analyzed in Ref ( Jonckheere et al\ 
2009b). The other parameters are eo and Vs = 0.2A. The 



insets show the comparison of the dc current and the corre- 
sponding harmonics for Fa^ = (black line) and Fa^ = 0.02 A 
(red line). Reprinted figure with permission from T. Jonck- 
heere e£_a/^hysical Review B 80, 184510, 2009 (J onckheere] 



Society. 



et al. 2009b). Copyright (2009) by the American Physical 




FIG. 52 Linear conductance for the side coupled dot be- 
tween a normal and a sup erconducting lead considered in 
Ref. (|Tanaka et al. 2008') as a function of the lateral dot 
level 62 and for different values of the corresponding Coulomb 
interaction parameter 1/2- The upper panel shows a schematic 
representation of the setup considered in this work. Reprinted 
figure with permissio n from Y. Tanaka et al, Physical Review 
B 78, 035444, 2008 ([Tanaka et alj [20081). Copyright (2008) 



by the American Physical Society. 



increase followed by a decrease. In this work it was also 
shown that the conductance in the electron-hole sym- 
metric case can be tuned to the unitary limit for fixed U 
by varying the coupling to the superconducting lead F^. 
This possibility is analogous to the one already discussed 
for the single SQDN system in Sect. |IV[ 



In a subsequent paper (Tanaka et a/., 2010|), the same 



authors considered the case of a double QD in series con- 
nected to a normal and a superconducting lead (see upper 
panel of Fig. 53) using the same theoretical approach. 



They first showed that for the case where the dot cou- 
pled to the normal lead is in the electron-hole symmetric 
condition the problem can be mapped into an effective 
normal two impurity Anderson model in terms of the 
Bogoliubov operators. Furthermore, this transformation 
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allows to calculate the conductance using the Friedel sum 
rule in terms of a phase-shift given by 



(43) 



where E = ^/e^ 



62 being the level for 



the dot coupled to the superconductor and Q = 
(< TiaTia > + < lL^2a >) , lia indicating the Bo- 
goliubov operators. The authors identified three different 
regimes. For small t and U < 2Ts a regime correspond- 
ing to the a local superconducting singlet is found while 
for U > 2Ts the Kondo singlet state is formed. On the 
other hand, for large t the antiferromagnetic coupling be- 
tween the dots dominates. The behavior of the conduc- 
tance in these three different regimes is summarized in 
Fig. 53 As can be observed, for sufficiently small U /F^ 
the conductance can reach the unitary limit, the maxi- 
mum displacing towards smaller t values as U increases 
and eventually for U /F^ ~> 2 the unitary limit cannot 
be reached. A surprising feature appears for U /F^ ~ 2 
where the conductance reaches the unitary limit for two 
different t values, indicated by the dashed rectangle in 
Fig.|53l 



3. Non-local Andreev transport through single or double 
quantum dots 

Double dots coupled to several normal and supercon- 
ducting leads are receiving recently considerable atten- 
tion in connection to the possibility of producing non- 
local entangled electrons from the splitting of Cooper 



pairs (Herrmann et al, 2010 Hofstetter et al. 2009 



Recher et al. , 2001D. The basic ide as were first put for- 



ward in Ref. (Recher et al., 2001) where the multiter- 



minal geometry of Fig. |54| was considered. In this con- 
figuration, when a Cooper pair is injected from the SC 
lead it can either be transmitted as a whole to one of 
the normal leads by a local Andreev process or split with 
each of the electrons in the pair transmitted to a differ- 
ent lead (corresponding to a non-local or crossed Andreev 
reflection process). The advantage of the DQD set up is 
twofold: on the one side it allows to tune independently 
the two dot levels and on the other hand Coulomb in- 
teractions could be used to favor the splitting processes 
compared to the local ones. While this issue is of a great 
current interest it goes beyond the scope of the present 
review. We would just mention several works address- 
ing the non-local Andreev transport involving quantum 
dots and multiterminal configura tions in Refs. (Eldridge 
et a/.l 120101 iFutterer et a/.l 120091). 



VII. CONCLUDING REMARKS 

In this review article we have summarized the most rel- 
evant published work related to superconducting trans- 
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FIG. 53 (a) Linear conductance for the series double quantum 
dot betw een a normal and a s uperconducting lead discussed 
in Ref. ( [Tanaka et al. |2010 | as a function of the interdot 
tunneling t and for different values of U /Vs. Panel (b) shows 
the corresponding induced pairing correlation in the dots re- 
gion. The results were obtained using the NRG method in 
the A ^ oo limit. Upper panel: setup considered in this 
work. Reprinted figure with permissio n from Y. Tanaka e t 
al, Physical Review B 81, 075404, 2010 ("Tanaka et aH|2010 ). 
Copyright (2010) by the American Physical Society. 



port in quantum dots systems. The large variety of topics 
that we have covered give an indication of the great ac- 
tivity which this field has shown in recent years. Due 
to the limited space it has become necessary to restrict 
somehow its content. For this purpose we have chosen to 
give priority to the more basic topics like Josephson ef- 
fect and Andreev transport through single level quantum 
dots, and had left aside some interesting but more spe- 
cialized situations like those involving ferromagnetic ma- 
terials or unconventional superconductors. In the same 
way, we have not covered in this review the response of 
these systems to exter nal ac fields, like phot on assisted 
transport in S-QD-S (Cho a/.[ [19991 |Zhu 2002) 



or N-QD-S systems ( |Zhao[ |1998[ fZhao and Wangl |2QQT| 



and adiabatic pumping in NDQS systems ( Splettstoesser 



et al. 2007). 



Even within the basic topics discussed in this review 
there remain several issues which are not completely un- 
derstood and deserve further analysis. Among them we 
may point out: 1) the confiicting description of ABs 
within the different approximation schemes for S-QD-S 
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als of combining these systems with topological insula- 



FIG. 54 Schematic representation of a double quantum dot 
coupled to a central superconducting and t wo lat eral normal 
leads in the proposal of Ref. (Recher et al. 2001 ) for produc- 



ing entangled electron pairs by means of non-local Andreev 
processes. Reprinted figure with permissi on from P. Recher e t 
al, Physical Review B 63, 165314, 2001 ([Recher et al 



Copyright (2001) by the American Physical Society 



2001) 



tors are still on an speculative level ( Golub et al. 



2011). 



We expect nevertheless that these issues would exhibit a 
great development in the near future. 
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